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ABSTRACT 


Worldview-2  spectral  imagery  acquired  over  Duck,  NC,  and  Pendleton,  CA,  were 
analyzed  to  extract  Bidirectional  Reflectance  Distribution  Functions  (BRDF)  for  8 
spectral  bands.  Spectral  data  in  the  visible  and  near- infrared  bands  were  acquired  for  15 
azimuth/elevation  values  during  the  Duck  NC  orbit  pass;  data  were  collected  at  ten- 
second  intervals.  Ten  images  were  acquired  over  Pendleton.  Orthoready  images  were 
coregistered  using  first-order  polynomials  for  the  two  image  sequences.  BRDF  profiles 
have  been  created  for  scene  elements:  vegetation,  asphalt,  sand,  and  water.  These  data 
allowed  for  unique  high-spatial  resolution  BRDF  profiles,  which  can  be  used  to  inform 
terrain  classification  and  target  detection  in  satellite  imagery. 
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I.  INTRODUCTION 


Classical  remote  sensing  systems  have  viewed  the  earth  from  a  perspective,  which 
has  been  limited  to  a  single  view,  generally  in  a  near-nadir  configuration.  Of  necessity, 
variations  in  reflectance  with  illumination  angle  and  view  angle  have  been  ignored. 
LANDSAT  is  a  classic  illustration  of  this  in  Earth  resource  systems.  Variations  in 
reflectance  with  angles  are  quantified  by  the  BiDirectional  Reflectance  Distribution 
Function,  abbreviated  “BRDF”.  This  quantity  has  been  measured  in  laboratory  settings, 
and  is  just  starting  to  be  observed  from  airborne  and  satellite  systems.  Satellites  offer  the 
best  way  to  make  regular  and  frequent  BRDF  measurements  on  a  large  scale.  With  a 
large  number  of  observations,  it  is  possible  to  accurately  model  the  global  BRDF.  This 
could  provide  benefits  to  applications  such  as  terrain  mapping.  This  thesis  addresses 
BRDF  observations  obtained  over  Duck,  NC  and  Pendleton,  CA,  utilizing  successive 
images  taken  by  the  Worldview-2  commercial  satellite  imaging  system.  The  advantages 
of  Worldview-2  over  its  predecessors  are  vastly  improved  spectral  and  spatial  resolution 
and  a  high  revisit  time,  which  serve  to  provide  more  data  points  to  the  model.  This  thesis 
displays  the  improvements  in  BRDF  observations  on  a  small-scale. 
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II.  BACKGROUND 


A.  REMOTE  SENSING 

Remote  sensing  involves  the  collection  of  data  from  sensors,  which  are  not  in 
direct  contact  with  the  object  being  observed.  Platforms  such  as  airplanes,  buoys,  and 
satellites  are  ideal  candidates  for  remote  sensing  instruments. 

Radiometry  is  the  field  concerned  with  the  measurement  of  the  electromagnetic 
spectrum.  Radiometry  deals  in  absolute  power,  unlike  the  field  of  photometry,  which 
measures  brightness/intensity  as  perceived  by  the  human  eye.  Radiometers  are  the  most 
common  remote  sensing  instrument  and  are  used  to  measure  the  radiant  flux,  or  power,  of 
electromagnetic  radiation.  Common  radiometry  terms  are  explained  in  the  next  section 
of  this  chapter. 

Many  products  are  derived  from  remote  sensing  instruments,  and  resolution 
properties  can  vary  widely  between  acquired  imagery.  Spatial,  spectral,  radiometric,  and 
temporal  resolution  are  limiting  features  of  the  specific  instruments  used  to  obtain 
imagery. 

Spatial  resolution  involves  a  sensor’s  ability  to  distinguish  physical  separation  of 
two  objects  viewed.  In  the  remote  sensing  community,  low  spatial  resolution  is  defined 
as  30  m  or  greater,  medium  ranges  from  4-30  m,  and  high  resolution  imagery  is  less  than 
4  m.  Excellent  satellite  imagery  is  currently  on  the  order  of  0.5  m  spatial  resolution. 


Figure  1.  Camp  Randall  Field,  Madison,  WI,  Spatial  Resolution  of  30m,  10m,  and 
lm  (Image:  University  of  Wisconsin,  Institute  for  Environmental  Studies) 
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Spectral  resolution  typically  refers  to  the  number  of  bands  over  the  spectral  range 
of  measured  frequencies.  For  example,  if  an  instrument  had  only  one  spectral  band  over 
the  visible  spectrum,  a  photon  of  red  light  (~700nm)  observed  by  the  instrument  would 
be  indistinguishable  from  a  measured  photon  of  green  light  (~500nm).  Panchromatic 
imagers  are  an  example  of  single-band  imagers.  However,  if  an  instrument  had 
numerous  incremental  bands,  then  the  instrument  could  distinguish  between  the  visible 
colors.  Multispectral  imaging  (MSI)  and  hyperspectral  imaging  (HSI)  refers  to  the  level 
of  spectral  resolution.  This  concept  is  illustrated  in  Figure  2. 


Figure  2.  Differences  between  Multi-spectral  imagery  and  Hyperspectral  imagery 


Radiometric  resolution  is  the  number  of  different  intensities  the  sensor  can  detect. 
Greater  resolution  allows  for  finer  distinguishability  between  intensities.  It  is  typically 
limited  by  the  available  number  of  bits.  For  example,  an  eight  digit  binary  string  might 
be  capable  of  defining  256  shades  of  gray;  however,  the  limiting  factor  on  radiometric 
resolution  is  often  noise  in  the  system. 
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The  final  type  of  resolution,  temporal  resolution,  is  a  description  of  how  often  a 
sensor  can  obtain  imagery  of  a  region  of  interest  (ROI),  also  known  as  its  revisit  time.  It 
is  useful  for  change  detection  in  imagery.  For  example,  deforestation  studies  do  not 
require  a  high  temporal  resolution,  as  change  occurs  over  a  long  period  of  time,  whereas 
the  intelligence  community  requires  high  temporal  resolution  as  well  as  low. 

B.  PHYSICS  OF  LIGHT 

While  the  term  Tight’  may  sometimes  be  used  interchangeably  with  ‘visible’, 
Tight’  is  a  general  term  for  all  wavelengths  and  frequencies  in  the  electromagnetic 
spectrum.  The  spectrum  is  illustrated  in  Figure  3.  The  visible  portion  of  the  spectrum 
refers  to  the  wavelengths  visible  to  the  human  eye,  approximately  400-700  mn.  All 
wavelengths/frequencies  possess  defining  characteristics  of  intensity,  polarization,  and 
phase,  and  where  necessary,  the  author  will  refer  to  wavelengths  by  their  general 
category  (i.e.,  visible,  near-infrared  (NIR),  ultraviolet,  etc.).  The  general  definitions  of 
these  categories  are  delineated  in  Table  1. 
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Figure  3.  The  Electromagnetic  Spectrum 
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Name 

Wavelength  (pm) 

Visible 

0.38-0.75 

Red 

0.62  -  0.75 

Orange 

0.59-0.62 

Yellow 

0.57-0.59 

Green 

0.495  -0.57 

Cyan 

0.476  -  0.495 

Blue 

0.45  -  0.475 

Violet 

0.38-0.45 

NIR 

0.75-1.4 

SWIR 

1.4 -3.0 

MWIR 

3.0 -8.0 

LWIR 

8.0-15.0 

FIR 

14.0-1000 

Table  1.  Visible  and  Infrared  Bands  defined  by  Wavelength 


1.  Interactions  With  Matter 

Light  can  be  considered  as  radiated  or  irradiated.  These  terms  refer  to  the  source 
and  destination  of  energy.  Radiation  refers  to  the  emission  of  energy  from  an  object, 
while  irradiation  refers  to  the  interaction  of  radiated  energy  upon  a  surface.  A  prime 
example  of  a  radiator  is  the  Sun,  while  the  Earth,  interacting  with  the  Sun’s  emitted 
energy,  would  be  considered  irradiated  by  the  energy.  Electromagnetic  radiation 
interacts  with  matter  in  one  of  four  ways:  scattering,  transmission,  absorption,  and 
reflection.  How  light  behaves  upon  striking  a  medium  depends  upon  the  characteristics 
of  the  medium,  the  angle  of  incidence,  and  the  characteristics  of  the  light  itself  [1]. 
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Figure  4.  Electromagnetic  Radiation  Interactions  with  Matter  [2] 


Scattering  does  not  change  the  energy  of  the  electromagnetic  wave  but  changes 
the  direction  of  the  wave.  The  interaction  causes  the  wave  to  be  redirected  in  a  myriad  of 
directions,  maintaining  wavelength,  but  with  each  resulting  wave  at  a  reduced  intensity  of 
the  original. 

Transmission  is  the  process  of  an  energy  wave  passing  through  a  medium.  It  is 
dependent  upon  the  properties  of  the  medium  the  energy  wave  is  passing  from,  properties 
of  the  medium  it  is  entering  into,  and  the  angle  of  incidence  between  the  two. 
Transmission  describes  how  a  wave’s  angle  of  incidence,  9i,  referenced  from  the  normal, 
is  affected  due  to  a  velocity  change  caused  by  its  passage  from  one  medium  to  another.  It 
results  in  a  velocity  and  direction  change  of  the  wave,  called  refraction.  This  concept  can 
be  illustrated  mathematically.  Each  medium  has  an  intrinsic  property,  called  an  index  of 
refraction,  r|,  which  is  described  in  Equation  (1): 
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(1) 


rj 


The  index  of  refraction  is  the  ratio  of  the  velocity  of  the  radiation  in  a  vacuum,  c,  to  the 
velocity  in  a  medium,  v.  As  the  energy  passes  from  one  medium  to  the  next,  the 
interaction  is  described  by  Snell’s  Law,  summarized  as  follows: 

sin  9X  _  Vj  _  n2 


sin  02  v2  /? 


(2) 


The  resulting  angle  of  refraction  is  02,  which  is  also  referenced  from  the  nonnal. 


Absorption  occurs  when  the  irradiated  object  does  not  allow  transmission  through 
the  material,  but  instead  the  energy  is  converted  into  internal  heat  that  is  then  reemitted  as 
thennal  radiation. 


Specular  reflection  involves  the  bounce  back  of  electromagnetic  radiation  by  a 
surface  upon  which  the  radiation  is  incident.  The  law  of  reflection  defines  a  plane,  the 
plane  of  incidence,  in  which  the  incident  wave,  reflected  wave,  and  surface  normal  lie. 
This  law  assumes  a  smooth  surface,  which  results  in  a  specular  reflection,  and  the  angle 
of  incidence  equals  the  angle  of  reflectance.  The  reflective  nature  of  a  surface  is 
summarized  by  the  equation: 


R  = 


Vj-Vj 


\2  f  A 


vVi+V2y 


(3) 


.rl2+rl,y 

Where  again  v  is  the  velocity  in  a  medium  and  q  is  the  index  of  refraction.  The 
reflectivity  of  the  air-water  interface  results  in  a  value  of  approximately  0.02,  or  2%. 
This  means,  2%  of  the  energy  is  reflected  back  and  98%  of  the  energy  is  transmitted 
through  the  water. 


2.  Solid  Angle 

Solid  angles,  quantified  by  the  definition  of  steradians,  are  useful  when  describing 
properties  and  behavior  of  light.  A  solid  angle,  f2,  refers  to  the  cone  of  angles  which 
subtend  an  area  on  the  surface  of  a  sphere,  the  area  of  which  is  equal  to  the  area  occupied 
by  an  object.  Mathematically,  solid  angles  are  defined  by  Equation  (4).  The  concept  is 
illustrated  by  Figure  5. 
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(4) 


While  a  radian  is  defined  in  2-D,  a  steradian  is  a  3-D  concept.  Steradian,  meaning 

solid  radian,  allows  solid  angles  to  be  quantified.  A  steradian  is  defined  on  a  sphere  of 

2  2 

radius,  R.  It  is  the  solid  angle  required  to  subtend  an  area  equal  to  R  (A=  R“).  When 
substituted  into  Equation  (4),  the  equation  evaluates  to  1 . 


Figure  6. 


Lengths  r 


Area  =  r2 


(a)  Radian  defined  (b)  Steradian  defined 


2 

The  surface  area  of  a  sphere  is  47iR  ,  regardless  of  the  size  of  the  sphere.  Utilizing 
Equation  (4),  the  solid  angle  of  any  sphere  is  calculated  as 
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(5) 


Calculating  steradians  removes  the  sphere’s  radius  from  consideration.  Steradians 
quantify  the  relative  amount  of  the  sphere  surrounding  a  point  that  is  occupied  by  an 
object  when  viewed  from  the  center  of  the  sphere.  An  object  that  obscures  the  entire 
region  surrounding  a  point  would  be  said  to  occupy  47t  steradians.  Conversely,  an  object 
that  appears  to  cover  only  half  of  the  viewable  sphere,  a  hemisphere,  occupies  2n  [sr] . 
Note  that  the  steradian  is  a  dimensionless  quantity,  and  [sr]  or  the  term  ‘steradian’,  are 
included  solely  for  clarity.  For  clarification,  when  an  object  subtends  only  a  portion  of  a 
given  sphere,  Equation  (4)  is  rewritten  as 

dn  =  d4  (6) 

R2 

3.  Foreshortening 

An  understanding  of  foreshortening  is  necessary  to  understand  the  subtle 
differences  in  radiometric  terms.  This  concept  is  clearly  explained  in  Figure  7. 


Edge  of  com 


A 


A  sees  the  diameter  of  the  coin  in  proportion 
to  the  length  of  (he  blue  line,  whereas  B  sees 
it  in  proportion  to  the  length  of  the  red  line. 


Bill  B  sees 


The  foreshortening  is  in  proportion  to  the 
cosine  of  the  angle  with  the  normal 
(next  slide)  n 


10 
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Figure  7.  Foreshortening  [3] 


Foreshortening  is  a  term  frequently  used  in  graphic  design,  and  it  describes  how 
an  image’s  appearance  varies  with  the  viewer’s  perspective. 

4.  Irradiance 

Irradiance  is  the  flux  per  unit  area  onto  a  surface  and  is  defined  as 


(V) 


where  <1>  is  the  radiant  flux,  or  power.  It  is  not  foreshortened.  Irradiance  implies 
directionality,  typically  power  inbound  to  a  surface  [4], 

5.  Radiance 

Radiant  intensity  describes  the  flux  per  unit  solid  angle  from  a  point  source  into  a 
particular  direction.  It  provides  directional  information  but  no  spatial  information,  and  it 
is  defined  as 


(8) 


where  df2  is  the  element  of  solid  angle.  Radiance,  L,  introduces  the  spatial  element  [4]. 
It  is  the  flux  per  unit  projected  area  per  unit  solid  angle.  As  it  deals  with  the  amount  of 
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energy  passing  through  a  specific  solid  angle  and  onto  an  area,  it  requires  the  introduction 
of  the  concept  of  foreshortening. 

L  = - — -  [Wsr'W  ]  (9) 

dCldA  cos  6 

Manipulating  Equations  (7), (8),  and  (9) 

dO 

L  = - — - =  — ^ —  =  — — —  (10) 

d Cld A  cos  0  cosdd Q  cos  OdLl 

Radiance  is  not  identical  across  all  wavelengths;  therefore,  a  tenn  is  needed  to  describe 
the  differing  quantity  with  respect  to  wavelength.  Introduction  of  the  term,  spectral 
radiance,  Lx,  allows  for  this  differentiation. 

6.  Bidirectional  Reflectance  Distribution  Function 

The  Bidirectional  reflectance  distribution  function  (BRDF)  describes  the 
scattering  characteristics  of  a  material  or  object  by  describing  how  light  from  a  source, 
incident  on  a  target,  is  reflected  in  a  given  direction.  The  BRDF  equation  is  the  ratio  of 
the  energy  scattered  by  the  target  in  a  particular  direction,  dependent  upon  wavelength, 
over  the  energy  that  was  incident  on  the  target  from  a  particular  direction  (BR)  integrated 
for  all  combinations  of  inbound  and  outbound  energy  directions  (DF).  BRDF  is 
dependent  on  both  the  angle  of  incidence  and  on  the  angle  of  reflectance.  These  angles 
are  defined  in  a  Cartesian  coordinate  system  by  a  polar  angle,  0,  measured  from  the 
surface  normal,  and  an  azimuthal  angle,  (p,  measured  from  the  x-axis.  The  geometry  is 
defined  in  Figure  8. 
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Figure  8.  Reflection  Geometry  used  in  the  definition  of  the  BRDF 


Symbol 

Name 

Unit 

Ei,Er 

Irradiance 

W 

m~ 

Ej  ,Ly 

Radiance 

W 

m2sr 

L\J,Lx,r 

Spectral  radiance 

W 

m3  sr 

f 

BRDF 

1 

sr 

Table  2.  Units  of  Radiometric  Terms 


The  BRDF  equation  is  [4] [5]  [6]: 


f/  (k,0i,</>i,0r,</>r) 


dEx,r(E.  Or.  (J)i-\  r  i  .. 

- - - —  [sr  1 

dEx,  i(k,  Oi,  (j)i) 


(ID 


Here,  dExj(k,Q j,(pi)  is  the  incident  spectral  irradiance  and  dLx,r(k,Qr,( pr)  is  the  reflected 
spectral  radiance  with  regards  to  wavelength.  This  equation  assumes  unifonn  irradiance 
over  a  large  area  of  a  uniform,  flat,  and  isotropic  surface.  Wavelength  is  disregarded  to 
simplify  the  equation.  Now  the  equation  is  written  in  terms  of  the  incident  irradiance  and 
reflected  radiance. 


£  (Oi,(f>i,Or,(f>r) 


dLr{Or,(j)r)  j 

dEi(Oi,</)i)  1  J 


(12) 
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Radiance,  unlike  irradiance,  does  not  involve  a  surface,  only  energy  at  a  point.  This 
allows  it  to  interchangeably  represent  energy  inbound  or  outbound.  To  elaborate, 
radiance  leaving  a  point,  p,  in  the  direction  of  a  point,  q,  is  the  same  as  radiance  arriving 
at  q  from  p.  This  concept  allows  the  BRDF  equation,  Equation  (11),  to  be  put  solely  in 
terms  of  radiance.  The  process  is  begun  by  manipulating  Equation  (10), 

dE  =  LcosddQ  (13) 

Then  Equation  (13)  is  rewritten  with  the  appropriate  reference  angles 

dEi(Qi,  (pi)  =  L  i{6i,  (pi)  cos  Od  Q,  ( 1 4) 

And  substituted  into  Equation  (12) 

fr  (  di,  (pi,  dr,  (pr)= - ^ -  [sf'  ]  (15) 

L/(&,  <pi)  cos  did  Eh 

Equation  (15),  another  way  of  expressing  the  BRDF  equation,  now  relates  solely  to 
radiance. 


7.  Anisotropy 

Anisotropy  is  the  property  of  being  directionally  dependent.  In  contrast,  an 
isotropic  surface  displays  a  property  identically  in  all  directions.  The  study  of  BRDF  is 
concerned  with  the  property  of  reflectance.  An  example  of  an  isotropic  reflector  is 
unfinished  wood,  as  it  appears  the  same  from  all  viewing  angles,  while  finished  wood 
displays  anisotropic  reflection.  The  anisotropic  nature  of  finished  wood  is  displayed  in 
Figure  9. 


Figure  9.  Anisotropic  nature  of  Finished  Wood 
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Most  land  surfaces  are  strongly  anisotropic  reflectors  in  the  optical  to  thermal 
wavelengths.  The  main  source  of  this  anisotropic  nature  is  the  three-dimensional 
character  of  a  surface.  Shadowing  and  image  obstruction  are  typical  reasons  an  object 
differs  in  its  appearance  when  viewed  from  different  angles.  This  manifests  in  changes  in 
surface  radiance.  An  accurate  BRDF  is  the  appropriate  method  for  capturing  and 
expressing  this  anisotropic  nature. 


Figure  10.  Anisotropy  of  Natural  Features.  Top:  Soybean  Field,  Bottom:  Black 
Spruce  Forest,  Left:  Backscattering  (sun  behind  observer),  Right: 
forwardscattering  (sun  opposite  observer),  Photograph  by  Don  Deering. 
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The  images  in  Figure  10.  ,  were  taken  as  part  of  the  Goddard  Space  Flight 
Center’s  PARABOLA  experiment,  and  display  how  the  sun  angle,  with  respect  to  the 
observer,  causes  a  difference  in  the  appearance  of  the  object,  illustrating  its  anisotropic 
nature.  Backscattering,  when  the  sun  is  behind  the  observer,  for  example,  causes  a  hot¬ 
spot  effect,  which  hides  all  shadows.  Forwardscattering,  where  the  sun  is  opposite  the 
observer,  shows  effects  such  as  specular  reflection  and  increased  shadowing  [7]. 

C.  GROUND-BASED  SENSORS 

Ground-based  sensors  are  required  to  measure  surface  BRDF  and  validate  the 
results  of  the  space-borne  sensors.  Gonioreflectometers  are  a  type  of  ground 
measurement  tools,  while  the  PARABOLA  is  a  unique  system  designed  with  the  mission 
of  surface  BRDF  measurements. 

A  gonioreflectometer,  as  the  name  implies,  was  designed  to  measure  the  angles  of 
incident  and  reflected  light.  It  provides  more  control  over  measurements  taken  in  situ  as 
it  consists  of  a  light  source,  the  material  to  be  observed,  and  a  sensor  to  measure  the 
reflected  light.  As  BRDF  is  a  four-dimensionally-dependent  concept,  the  original 
gonioreflectometer,  proposed  in  Murray-Coleman  and  Smith’s  1990  article  for  the 
Journal  of  the  Illuminating  Engineering  Society,  provided  four  degrees  of  mechanical 
freedom.  Gonioreflectometers  have  provided  enormous  amounts  of  data  as  they  present  a 
way  to  accurately,  though  not  expeditiously,  measure  the  BRDF  over  a  large  combination 
of  angles  [8]. 

The  first  major  renovation  to  the  original  design  was  proposed  by  the  Lawrence 
Berkeley  Laboratory  in  1992.  It  differed  from  the  original  design  as  it  incorporated  a 
hemispherical  mirror  and  a  fish-eye  lens,  which  allowed  for  simultaneous  measurement 
of  all  angles  of  reflected  energy  at  a  given  incident  angle  without  movement  of  the 
camera.  Though  a  dramatic  simplification  of  the  BRDF  measurement  collection  process, 
the  design  was  not  without  its  own  flaws,  such  as  distortion  caused  by  the  fish  eye  lens 
and  unintended  re-illumination  of  the  sample. 
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Source  Dr  v^r 


Figure  11.  Murray-Coleman  and  Smith’s  Gonioreflectometer  (1990) 

Side  Vita 


Figure  12.  Lawrence  Berkeley  Laboratory’s  Gonioreflectometer  (1992)  [9] 
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Figure  13.  Incidence-plane  BRDF  of  rough  aluminum  surface  for  several  incidence 

angles,  0;  [10] 

Figure  13.  displays  a  sample  of  measurements  taken  from  the  results  of  a  study 
into  automation  of  a  three-axis  gonioreflectometer.  The  figure  is  taken  from  a  full 
hemispherical  dataset.  It  displays  a  rise  in  the  BRDF  value  when  the  incident  angle  and 
the  reflected  angle  coincide  [10]. 

PARABOLA,  the  Portable  Apparatus  or  Rapid  Acquisitions  of  Bidirectional 
Observations  of  Land  and  Atmosphere,  is  a  three  channel,  rotating  head  radiometer.  It 
has  three  detectors,  two  silicon  and  one  germanium  solid  state,  with  filters  corresponding 
to  the  spectral  bands,  630-690,  760-900,  and  1550-1750  mn,  respectively.  The  two-axis, 
two-motor  rotation  of  the  head  enables  a  near-complete  sampling  of  the  entire  sky/ground 
sphere  with  a  FOV  of  15°.  During  its  3-minute  cycle  it  measures  72  azimuthal  and  37 
elevation  positions,  generating  radiance  measurements  for  both  the  ground  and  sky.  Its 
primary  purpose  is  to  research  the  nature  of  BRDF  and  was  designed  to  provide  in-field 
verification  of  MISR  data. 
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Figure  14. 


Portable  Apparatus  or  Rapid  Acquisitions  of  Bidirectional  Observations 
of  Land  and  Atmosphere  (PARABOLA)  Ground  Sensor  [11] 


□  PARABOLA  - 4-scale  model 


View  Zenith  Angle 


Figure  15.  Reflectance  values  for  multiple  solar  zenith  angles  in  Red  (Top)  and 
NIR  (Bottom)  bands  of  PARABOLA  and  Chen-Leblanc’s  radiative  transfer 

model  are  compared 
19 


Experiments  were  done  in  the  forests  in  Saskatchewan,  Canada,  utilizing 
PARABOLA.  PARABOLA  data  was  compared  to  results  calculated  by  a  geometric 
optical  model  radiative-transfer  model.  As  Ligure  15.  shows,  the  model  tested  simulates 
the  measurements  closely,  with  the  exception  of  at  large  view  zenith  angles  on  the 
backscattering  side  [12]. 

D.  SPACE-BORNE  SENSORS:  MULTI-ANGLE  AND  MULTIPLE  ANGLE 

SATELLITES 

Satellites  have  long  been  used  to  make  measurements  of  the  Earth.  The  NASA 
concept  called  the  Earth  Observing  System  (EOS)  consists  of  a  series  of  coordinated 
polar  and  low-inclination  satellites  designed  for  global  observations.  Most  of  these 
satellites  are  not  in  a  position  to  make  BRDF  measurements,  as  they  generally  offer  only 
a  nadir  view.  A  small  number  of  sensors  have  begun  to  offer  routine  multiple  angle 
views.  The  advent  of  IKONOS  in  1999  allowed  for  significant  off-nadir  collection.  This 
raises  issues  on  how  to  deal  with  off-nadir  views  in  the  analysis  of  spectral  data.  These 
satellites  aid  in  atmospheric  corrections  as  they  make  timely  and  successive,  and  in  the 
case  of  multi-angle,  simultaneous  measurements  of  the  atmospheric  aerosols  which 
informs  and  improves  atmospheric  modeling  [13].  The  use  of  cross-track  and  along-track 
instruments  counters  the  narrow  field  of  view  (FOV)  argument  and  allows  for  an  increase 
in  the  range  of  angular  observations  that  can  be  made. 

Multi-angle  satellites  acquire  multiple  images  of  the  same  object  at  differing 
angles  through  the  simultaneous  employment  of  multiple  sensors.  These  systems  allow 
for  instantaneous  sampling  of  BRDF.  Multiple  angle  satellites  acquire  multiple  images 
of  the  same  object  at  differing  angles  through  use  of  a  single  sensor  over  a  period  of  time. 
MISR  and  POLDER  are  examples  of  multi-angle  systems,  while  PROBA/CHRIS  is  a 
multiple  angle  imager.  MISR  was  originally  launched  aboard  the  Terra  (EOS  AM) 
satellite.  Images  from  multi-angle  and  multiple  angle  satellites  are  advantageous  to  the 
study  of  BRDF  as  they  can  be  used  to  validate  models  or  be  coupled  with  models  to 
provide  a  more  accurate  picture  of  BRDF. 
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1. 


Multi-angle  Imaging  SpectroRadiometer  (MISR) 


The  Jet  Propulsion  Laboratory  (JPL)  designed  the  MISR  sensor,  which  was 
launched  on  NASA’s  EOS  flagship  satellite,  Terra,  in  December  of  1999.  MISR  has  nine 
cameras  that  simultaneously  image  the  Earth  at  nine  discrete  angles.  Cameras  point 
forward  and  aftward  of  the  local  vertical  at  26.1°,  45.6°,  60.0°,  and  70.5°,  and  at  Nadir. 
MISR  images  in  the  blue,  green,  red,  and  near-infrared  spectral  bands,  and  its  spatial 
resolution  ranges  from  1  km  in  its  normal  operating  mode  to  275  m  in  local  mode.  MISR 
is  a  pushbroom  sensor,  also  called  a  cross-track  imager,  which  means  its 
image  is  obtained  through  the  motion  of  the  satellite  itself.  The  counter  type,  a 
whiskbroom  sensor  or  along-track  imager,  acquires  its  images  through  the  additional 
side-to-side  sweeping  motion  of  the  sensor  [e.g.,  MODIS]. 

The  multi-angle  capability  of  MISR  aids  in  furthering  understanding  of  how 
sunlight  is  scattered  under  varying  conditions.  Also,  it  provides  information  regarding 
monthly,  seasonal,  and  long-term  trends  different  types  of  clouds,  aerosol  particles,  and 
surfaces.  More  specifically,  it  allowed  for  advanced  study  of  the  amount,  type,  and 
source  of  aerosol  particles;  characteristics  of  clouds;  and  land  cover  distribution. 


Height:  3.5  meters 

Length:  6.8  meters 

Weight  5,1 90  kilograms 

Power:  2,530  watts  (average) 

Instrument  Data  Rate:  1 8,545  kilobytes 

per  second  (average) 

Design  Lifetime:  6  years 


CERES 


Figure  16.  Terra  Spacecraft  with  5  Climate-Monitoring  Sensors:  MODIS,  MISR, 

ASTER,  MOPITT,  and  CERES 
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D  view  angles  at  Earth 
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Calibrated  measurements 
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Backward -vlawl  rig 
cameras 


Figure  17.  Multi-angle  Imaging  SpectroRadiometer  (MISR)  observes  the  Earth 
simultaneously  with  nine  discrete  cameras  from  ±70.5°  in  four  spectral 

bands [14] 


Table  3.  MISR  Spectral  Bands 


Name 

Center  Wavelength  (mn) 

Blue 

446 

Green 

558 

Red 

672 

NIR 

867 

Numerous  studies  have  been  done  into  the  use  of  MISR  data  in  the  mapping  of 
BRDF.  The  improved  spatial  resolution  in  local  mode  and  the  multi-angle  viewing 
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capabilities  have  proven  beneficial.  In  a  study  over  the  Queensland,  Australia,  the  non¬ 
linear  Rahman-Pinty-Verstraete  (RPV)  model  and  the  linear  Ross-Thick  Li-Sparse 
Reciprocal  (Ross-Li)  model  were  compared  using  local  mode  MISR  data.  These  models 
are  further  explained  in  the  next  section. 


0.3 


Blue  Green  Red  NIR 


Li  RPV  near  hotspot  Hi  Ross-Li  near  hotspot 

□  RPV  mean  □  Ross-Li  mean 

Figure  18.  Mean  and  hot  spot  root  mean  squared  error  (RMSE)  statistics  displaying 

accuracy  of  RPV  and  Ross-Li  model  inversions 

Finding  the  RPV  model  a  better  fit,  it  was  compared  to  the  foliage  projective 
cover  (FPC)  and  it  was  assessed  that  the  resulting  BRDF  model  were  sufficient  to  map 
the  vegetation  in  the  grassland  and  woodland  region  of  the  Southern  Brigalow  Belt.  The 
study  concludes  that  the  use  of  Local  MISR  data  removes  the  need  to  use  complex 
atmospheric  radiative  transfer  algorithms  (RTAs)  within  the  scope  of  the  limited 
geographic  region  [15].  It  was  proposed  that  MISR  data  could  be  used  to  improve  the 
BRDF  and  albedo  products  created  by  the  Moderate  Resolution  Imaging 
SpectroRadiometer  (MODIS)  instrument,  also  onboard  the  Terra  satellite.  Preliminary 
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work  showed  a  reduction  in  bidirectional  reflectance  factors  (BRF)  by  10%  in  the  red  and 
green,  when  compared  to  MODIS-only  data  [16]. 

2.  Polarization  and  Directionality  of  the  Earth’s  Reflectances 
(POLDER) 

POLDER  was  launched  aboard  the  sun-synchronous,  polar-orbiting  Advanced 
Earth  Observing  Satellite  (ADEOS)  in  August  of  1996  and  POLDER  II  was  launched  on 
ADEOS  II  in  December  of  2002,  which  flew  until  2003.  ADEOS  and  its  instruments 
were  also  part  of  NASA’s  EOS  mission.  The  POLDER  instruments  were  important  to 
the  study  of  BRDF  measurements  as  they  provided  the  first  opportunity  to  sample  the 
BRDF  of  every  point  on  Earth  for  viewing  zenith  angles  up  to  60°  for  the  full  azimuth 
range,  at  a  spatial  resolution  of  approximately  6  km.  This  was  accomplished  through  the 
use  of  a  wide  FOV  lens  and  a  rotating  filter  and  polarizing  wheel  [17]  [18]. 

The  lens  allowed  for  a  maximum  FOV  of  114°,  specifically  ±43°  along  track  and 
±5 1  °  across  track.  The  polarizing  wheel  rotated  and  scanned  eight  narrow  spectral  bands 
in  the  visible  and  NIR  in  4.9  seconds.  The  surface  target  was  viewed  up  to  14  times 
during  the  satellite  overpass  with  various  viewing  angles.  Over  successive  days,  with  the 
assumption  of  cloud-free  days,  orbital  shifts  allowed  a  sampling  of  the  BRDF  over  the 
whole  instrument  field  of  view. 
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Figure  19.  Illustration  of  the  POLDER’s  rotating  polarizing  wheel.  Scanned  bands 
are  in  the  VNIR  bands  (A,=443,  490,  564,  670,  763,  765,  865,  and  910  nm). 

POLDER  filled  the  polarized  reflectance  observation  need.  An  additional  benefit 
of  POLDER  and  ADEOS  was  that  they  did  not  operate  solely  cross-track,  as  previous 
satellites  had.  This  allowed  for  less  restrictive  BRDF  measurements.  The  processing  of 
POLDER’s  data  was  improved  as  well.  POLDER’s  scientists  used  a  corrected  Rayleigh 
scattering  and  an  improved  atmospheric  scattering  model,  which  no  longer  relied  on  a 
simple  Lambertian  surface.  Analysis  of  the  data  in  follow-on  reports  ( [  1 9]  [20]  [2 1  ]), 
showed  that  the  imaging  systems  captured  the  hot  spot  and  specular  features  well. 
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Figure  20.  Sample  BRDF  product  derived  from  the  space-bome  POLDER 

instrument.  The  angular  space,  annotated  (p,  represents  the  azimuthal  delta 
between  the  sun  and  the  viewing  direction.  BRDF  is  at  a  maximum  in  the 
backscattering  direction  (hotspot,  cp=0°)  and  at  a  minimum  in  the  forward 
scattering  direction  (cp=180°)  [22] 


As  shown  in  Figure  20.  ,  BRDF  is  plotted  for  two  wavelengths,  670  nm  and  867 
nm.  The  BRDF  peaks  when  the  sun  and  viewing  directions  coincide.  This  is  evident  in 
the  principle  plane  displayed  in  Figure  21. 
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Figure  21 .  Reflectance  values  for  sample  region,  Tree  cover,  broadleaved, 
evergreen,  acquired  by  space-bome  POLDER  instrument  [22] 


3.  Compact  High  Resolution  Imaging  Spectrometer  (CHRIS) 

The  European  Space  Agency  (ESA)  launched  CHRIS  aboard  the  ‘Project  for  On- 
Board  Autonomy’  (PROBA-1)  satellite  in  2001.  The  instrument  provides  large  images, 
13  km  x  13  km,  at  medium  resolution,  17  m,  in  13  narrow  spectral  bands.  When 
resolution  is  reduced  to  34  m,  it  is  capable  of  imaging  in  62  spectral  bands.  The  sensor  is 
deemed  ‘multiple  angle’  because  the  satellite  has  four  reaction  wheels.  These  reaction 
wheels  allow  PROBA-1  to  slew  along-track,  which  gives  CHRIS  the  ability  to  image  a 
site  five  times  during  a  single  pass,  and  across-track,  which  reduces  its  revisit  time  of  a 
ROI  to  less  than  a  week.  The  high  spatial  resolution  of  CHRIS  allowed  for  validation  of 
other  space-bome  sensors  such  as  POLDER,  MODIS,  and  MISR. 
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Figure  22.  Compact  High  Resolution  Imaging  Spectrometer  (CHRIS)  [23] 

Research  has  suggested  the  usefulness  of  CHRIS  in  deriving  leaf  area  indices 
from  BRDF  model  inversion  [24].  Additional  research  into  the  usefulness  of  automatic 
land  cover  mapping  products  generated  from  CHRIS-derived  BRDF  data  utilizing  neural 
network  algorithms  has  proven  promising  [25],  and  improvements  to  fitting  the  retrieved 
CHRIS  data  have  been  proposed  through  the  use  of  an  atmospheric  correction  algorithm 
specific  for  PROB A/CHRIS  over  land  [26]. 

Figure  23.  displays  reflectance  data  acquired  by  PROBA/CHRIS.  It  displays  the 
reflectance  characteristics  of  four  different  regions  over  five  different  angles,  with 
vegetated  areas  depicted  in  the  left  two  charts,  and  man-made  materials  in  the  two  charts 
on  the  right.  The  vegetation  IR  ledge  can  be  clearly  observed  over  all  angles  in  the 
cornfield.  The  four  charts  illustrate  the  varying  dependence  of  the  materials  on  the  fly-by 
zenith  angle  (FZA),  the  viewing  angle.  The  curves  are  also  flatter  across  the  man-made 
materials,  indicating  a  lesser  dependence  on  wavelength.  The  multiple  angles  plotted 
reveal  the  dependence  of  the  material  on  viewing  angle.  Each  of  the  four  regions 
displays  different  angle-dependency. 


28 


Urban  fabrics 


U 

u 

c 

□ 


o 

a 


o 

CL 

o 


0.5 

0.4 

0.3 

0.2 

0.1 

0.0 

500  600  700  800  900  1000 

Wavelength  (nm) 


Industrial  units 


Figure  23.  Reflectance  data  derived  from  CHRIS/PRO  BA  over  a  range  of  zenith 

angles  [25] 


E.  BRDF  MODELS 

Given  the  complex  nature  of  measuring  BRDF,  approximations  have  been  made 
in  the  past  through  the  use  of  models.  These  models  have  been  extensively  developed  for 
use  in  3-D  graphics  and  modeling  [27].  These  models  are  also  tools  used  by  scientists  to 
improve  the  accuracy  of  calculations  and  to  provide  insight  into  affects  seen  in  acquired 
imagery.  These  models  are  constantly  being  updated  as  the  communities  understanding 
of  BRDF  increases.  Newer  models  are  inherently  more  complex.  Some  attempt  to 
account  for  more  variables,  while  others  focus  on  better  modeling  smaller  regions  of 
interest.  The  Lambertian  model  was  among  the  first  BRDF  models.  Torrance-Sparrow 
and  Oren-Nayar  are  other  follow-on  models  that  attempt  to  better  model  surface  BRDFs. 
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1. 


Lambertian  Model 


The  Lambertian  model  is  one  of  the  earliest  and  simplest  models.  It  was  derived 
by  Johann  Lambert  in  1760.  The  Lambertian  model  assumes  that  the  surface  is  perfectly 
isotropic  and  therefore  the  apparent  radiance  is  the  same  when  viewed  from  any  angle. 
Rough  surfaces  are  often  approximated  as  Lambertian  surfaces.  Lambert’s  cosine  law, 
also  known  as  the  cosine  emission  law  or  Lambert’s  emission  law,  is  the  mathematical 
explanation  for  this  postulation.  It  assumes  that  the  radiant  intensity,  I,  from  an  isotropic 
reflector  is  directly  proportional  to  the  cosine  of  the  angle  between  the  observer’s  line  of 
sight  and  the  surface  nonnal.  The  radiance  is  constant  because  as  the  intensity  is  reduced 
as  the  distance  increases,  so  too  does  the  observed  area.  See  Figure  24.  This  model  is 
for  a  perfectly  diffused  object.  It  is  often  used  to  represent  the  diffuse  element  in  newer, 
more  complex  models.  The  BRDF  of  a  Lambertian  surface  is  a  constant. 


Figure  24.  Explanation  of  Lambert’s  Cosine  Law.  Io  represents  the  observed 

radiance. 


In  an  experiment  done  on  South  African  savanna  vegetation,  researchers  found 
that  a  leaf,  when  intact  displayed  specular  reflection,  but  when  crushed  more  closely 
resembled  a  diffuse  surface.  As  the  particle  size  decreased  with  further  crushing,  the 
BRDF  of  the  leaf  flattened  toward  a  constant,  displaying  a  BRDF  of  a  near-Lambertian 
surface.  This  is  displayed  in  Figure  25. 
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(a) 


Viewing  Angle  (°) 


(b) 


Viewing  Angle  (°) 


Figure  25.  Laboratory-based  BRDF  of  whole,  4-,  and  2 -mm  cut  poplar  leaves  at  (a) 
normal  incidence  and  340  mn  and  (b)  67°  incidence  and  860  mn  [28] 

2.  Rahman-Pinty-Verstraete  Model 

The  semi-empirical  RPV  model  is  dependent  on  three  terms,  rhoO,  k,  and  0, 
described  by  the  function: 

p  s(0i,<pi\0r,(pr,X)  =  rhoO*M(#;,  dr,  k)*F(g)*[  1+R(G)]  (16) 

where  rhoO,  the  intensity  of  the  surface  reflectance,  and  k,  the  level  of  anisotropy  of  the 
surface,  are  empirical  surface  parameters,  and  0  controls  the  function  F(g)  and 
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influences  the  relative  forward  and  backscattering.  The  model  expresses  the  reflectance 
ps,  in  terms  of  illumination  from  a  direction,  0i  and  (pi,  and  the  direction  of  observation, 
0r  and  cpr  [29]. 


3.  Ross-Thick  Li-Sparse  Model 

The  Ross-Li  model  is  a  linear,  kernel-driven  model  where  BRDF  is  expressed  as  a 
sum  of  theoretically  constructed  kernel  functions,  dependent  upon  solar  zenith  angle,  0S, 
viewing  zenith  angle,  0V,  and  relative  azimuth  angle,  cp,  between  the  sun  and  viewer 
planes.  The  equation  describing  the  model  is: 

p  {9  s,  Ov,  (p)  =  ao  +  ai*  f  i  (6f ,  Ov,  <p)  +  ai*  f  i  (Os,  Ov,  cp)  (17) 

The  two  kernels  that  make  up  the  sum  represent  single  scattering  from  a  dense 
leaf  canopy  (fi)  and  sparse  surface  objects  with  geometric-optical  mutual  shadowing  (f2). 
Also  included  are  coefficients  that  represent  isotropic,  volumetric  and  geometric 
reflectance  a0,  &\,  and  a2  [30]. 

4.  Torrance-Sparrow  Model 

The  Torrance-Sparrow  model  was  created  in  1967.  It  was  one  of  the  first  models 
to  capture  off-specular  peak.  It  also  introduced  the  geometric  attenuation  factor,  G, 
which  allowed  for  masking  and  shadowing  [31][8]. 

As  the  model  was  based  in  geometrical  optics,  the  physics  of  light’s  motion  was 
simplified  to  consider  light  moving  as  a  ray,  vice  the  true  wave  nature  of  light.  This 
assumption  requires  the  size  of  the  wavelength  to  be  significantly  smaller  than  the  matter 
it  is  interacting  with.  The  model  surface  was  made  up  of  many  microfaceted  surfaces, 
considered  to  be  perfect  specular  reflectors,  with  differing  normals.  The  normals 
followed  a  Gaussian  probability  distribution.  The  model  assumed  the  diffuse  element 
was  due  to  internal  scattering  and  multiple  microfacet  reflections.  The  format  of  the 
model’s  specular  element  was  derived  using  Fresnel  reflectance  and  it  was  summed  with 
the  diffuse  element  to  achieve  the  end  equation. 
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The  model  fit  the  data  well,  however,  to  do  so,  it  was  necessary  to  input  several 
constants.  Its  simplified  physics  and  the  necessity  of  these  constants  are  the  major 
shortcomings  of  the  model. 


(c)  I  Jiwrcncc 
<PSNK=30.3722> 


<fl  Ward 

<PSNK=22.25W> 


(g)  Ward-Duer  <h)  Our  Anisotropic 

(PSNR=31.71 16)  •  PSNR=32. 1 210) 


Figure  26. 


Various  BRDF  models  plotted  against  a  Reference  Image  [32] 


A  follow  on  model,  the  Cook-Torrance  model  was  developed  from  the  Torrance- 
Sparrow  model,  and  offered  improvement  through  the  accountability  of  wavelength  and 
color  shifting.  Various  models  are  plotted  against  a  reference  image  in  Figure  26.  The 
Cook-Torrance  model  is  depicted  in  part  (c)  of  the  figure.  The  intensity  of  the  insets 
represents  disparity  between  the  model  and  the  reference  image. 

5.  Oren-Nayar  Model 

This  model  is  similar  to  the  Torrance-Sparrow  model  in  that  it  uses  a 
microfaceted  surface;  however,  it  differs  because  it  assumes  these  microfacets  are 
perfectly  diffuse  vice  perfectly  specular.  The  Oren-Nayar  model  was  created  in  1993, 
and  its  main  advantage  was  an  improvement  over  the  Lambertian  technique  for  modeling 

rough,  diffuse  surfaces.  Lambert’s  model  does  not  account  for  surface  roughness.  Rough 
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surfaces  have  varying  nonnals  across  the  surface  and  therefore  foreshortening  occurs 
disproportionally  between  the  intensity  and  area.  This  means  the  radiance  is  not  held 
constant  as  it  is  in  the  Lambertian  model  [33]. 


Real  Image  Lambertian  Model  Oren-Nayar  Model 


Figure  27.  Oren-Nayar  Model  for  Rough  Surfaces.  Note  the  improved  accuracy  of 

the  model  over  that  of  the  Lambertian  model. 
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III.  EXPERIMENTAL  APPROACH 


A.  PROBLEM  DEFINITION 

An  area  that  has  had  limited  research  attention  involves  the  accurate  calculation 
and  modeling  of  bidirectional  reflectance.  The  difficult  and  specific  nature  of  the 
calculations  has  driven  researchers  to  labor  to  improve  models  to  fit  the  limited  amount 
of  data.  The  calculation  of  BRDF  is  extremely  difficult  because  it  requires  individual 
measurements  of  all  combinations  of  direction  of  incoming  irradiance  from  a  source  and 
the  subsequent  outgoing  radiance  from  a  target  over  a  broad  spectrum  of  wavelengths. 
Measurements  have  been  done  in  the  field,  which  have  provided  precise  information 
regarding  a  few  specific  surfaces.  From  such  measurements,  models  have  been  created, 
which  may  provide  a  workable  representation  of  the  BRDF  under  similar  conditions  of 
material  and  terrain.  Tradeoffs  are  made  in  accuracy  in  an  effort  to  simplify  analysis,  as 
calculating  BRDF  for  each  experiment  would  be  tedious  and  time  consuming.  This  is  an 
acceptable  solution  for  some  but  a  poor  solution  for  others.  Further  work  is  needed  in 
finding  a  expeditious  and  cost-effect  method  to  retrieve  the  bidirectional  reflectance  over 
a  broad  spectrum  of  wavelengths,  angles,  and  geography. 

The  goal  of  this  thesis  is  to  illustrate  the  use  of  Worldview-2,  specifically  the 
advantage  provided  by  its  contribution  to  the  library  of  BRDF  data,  through  its  generation 
of  high  spatial  resolution  data  in  eight  wavelengths. 

B.  EQUIPMENT 

1.  Worldview-2 

The  data  for  analysis  was  collected  by  the  Worldview-2  satellite.  The 
Worldview-2  satellite  was  launched  by  DigitalGlobe  in  2009.  The  addition  of  this 
satellite  added  unique  remote  sensing  capabilities  that  were  previously  unavailable. 
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Figure  28.  Artist  rendering  of  Worldview-2  Satellite 


Four  spectral  bands  are  typically  provided  in  multispectral  imagery.  These 
include  bands  in  the  red,  green,  blue  and  near-infrared  spectrum.  These  bands  are 
provided  by  Worldview-2  and  by  its  predecessor,  Quickbird,  which  is  also  a  multispectral 
imager.  Worldview-2  differs  because  it  is  capable  of  imaging  in  four  new  bands.  A 
visual  and  numerical  definition  of  these  bands  is  provided  in  Figure  29.  and  Table  4.  , 
respectively.  The  first  new  band,  ‘coastal  blue’,  with  wavelength  limits  of  400-450  nm, 
offers  improved  bathymetry,  vegetation,  and  atmospheric  correction  study.  The  second 
additional  band  is  a  yellow  band,  with  a  range  of  585-625  nm,  which  aids  in  creating 
true-color  images  and  improves  target  characteristics  important  in  vegetation  analysis. 
The  third  band  is  a  red-edge  band.  Imaging  from  705-745  nm,  this  band  is  useful  for 
evaluating  the  health  of  vegetation  and  the  production  of  chlorophyll.  Finally, 
Worldview-2  provides  a  second  near- infrared  band,  capable  of  viewing  wavelengths  from 
860-1040  nm.  This  band  provides  an  overlap  in  the  original  near-infrared  band,  and  is 
less  susceptible  to  atmospheric  effects. 
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Relative  Response 


The  8  spectral  bands  of  WorldView-2 


Pan 

Coastal 

Blue 

Green 

Yellow 
Red 
Red  Edge 
NIRi 
NIRj 


3  SO  4S0  5  SO  650  750  BSO  950  1050 

Wavelength  (nm) 


WV02  Relative  Spectral  Radiance  Response 


Figure  29. 


Worldview-2  Bands  and  Relative  Spectral  Radiance  Response 
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Table  4.  The  Spectral  Bands  of  Worldview-2  [34] 


Band 

Center  (rnn) 

Min  (nm) 

Max  (nm) 

Name 

PAN 

632 

464 

801 

Panchromatic 

1 

427 

400 

450 

Coastal  Blue 

2 

478 

450 

510 

Blue 

3 

546 

510 

580 

Green 

4 

608 

585 

625 

Yellow 

5 

659 

630 

690 

Red 

6 

724 

705 

745 

Red-Edge 

7 

831 

770 

895 

NIR1 

8 

908 

860 

1040 

NIR2 

An  additional  benefit  of  Worldview-2  over  DigitalGlobe’s  previous  spectral 
imagery  satellite,  Quickbird,  is  improved  spatial  resolution.  Worldview-2  provides  a 
nominal  0.5  m  panchromatic  and  2.0  m  multispectral  imagery. 


Satellite 

Panchromatic  GSD  (m) 

Multispectral  GSD  (m) 

Nadir 

Off-Nadir 

Nadir 

Off-Nadir 

Quickbird 

0.61 

0.72  (25°) 

2.44 

2.88  (25°) 

Worldview- 1 

0.50 

0.55  (20°) 

N/A 

N/A 

Worldview -2 

0.46 

0.52  (20°) 

1.80 

2.40  (20°) 

Table  5.  DigitalGlobe  Satellites’  Panchromatic  and  Multispectral  Ground  Sample 

Distance  (GSD)  Resolution 


2.  The  Environment  for  Visualizing  Images  4.7  (ENVI®) 

A  number  of  different  tools  are  available  to  process  and  analyze  geospatial 
imagery.  Image  processing  done  in  support  of  this  thesis  was  done  using  ITT  Visual 
Infonnation  Solution’s  ENVI  software.  The  integrated  Interactive  Data  Language  (IDL) 
software  was  utilized  to  process  and  analyze  the  25  multispectral  image  products 
received. 

ENVI  provides  preprocessing  software  for  Worldview-2  imagery.  This  software 
allows  conversion  of  digital  numbers,  which  relate  relative  brightness  throughout  the 
image,  into  radiance.  This  is  done  with  a  known  gain  and  offset  correction,  embedded 
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within  the  imagery  files.  Alignment  and  spatial  correction  of  the  image  is  done  through 
use  of  ENVI’s  mosaicing  and  registration  tools.  Analysis  tools,  such  as  the  Region  of 
Interest  tool  and  the  Spatial  Profile  tool  were  also  used  to  extract  information  from  the 
imagery.  Further  explanation  of  the  employment  of  these  tools  is  provided  in  the 
Methodology  chapter. 
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IV.  METHODOLOGY 


A.  PRODUCT 

Multi-angle  data  were  acquired  as  part  of  an  experiment  in  bathymetry  derivation 
obtained  from  observing  wave  motion  [35].  Sites  were  chosen  for  advantageous  access 
and  topography.  The  utility  for  our  purposes  was  serendipitous.  Two  locations  were 
imaged:  Pendleton,  CA,  and  Duck,  NC.  The  latter  site  is  advantageous  for  our  work 
because  of  very  low  relief,  simplifying  registration. 

Imaging  by  Worldview-2  was  done  over  both  locations  on  March  24,  2010.  Ten 
images  were  taken  over  Pendleton,  CA,  and  15  over  Duck,  NC,  by  both  the  panchromatic 
imager  and  the  multispectral  imager.  The  multispectral  data  is  the  focus  of  analysis  in 
this  thesis. 


41 


Name 

Acquisition  Time 
(GMT) 

Mean  Satellite 
Azimuth  (deg) 

Mean  Satellite 
Elevation  (deg) 

1/P006 

15:41:35 

57.9 

46.9 

2/P004 

15:41:46 

63.4 

49.4 

3/P005 

15:41:57 

69.8 

51.6 

4/P013 

15:42:08 

77 

53.5 

5/P009 

15:42:18 

85 

54.9 

6/P003 

15:42:29 

93.8 

55.8 

7/P008 

15:42:40 

102.9 

56 

8/P010 

15:42:51 

112 

55.4 

9/P012 

15:43:01 

120.4 

54.3 

10/P014 

15:43:12 

128.1 

52.6 

1 1/P00 1 

15:43:22 

134.8 

50.6 

12/P015 

15:43:33 

140.5 

48.3 

13/P011 

15:43:43 

145.3 

46 

14/P002 

15:43:53 

149.5 

43.7 

15/P007 

15:44:03 

153 

41.3 

Name 

Solar 

Azimuth 

(deg) 

Solar 

Elevation 

(deg) 

Solar  Zenith 
Angle  (deg) 

Separation 
Angle  (deg) 

1/P006 

144.5 

50.1 

39.9 

66.61301 

2/P004 

144.5 

50.1 

39.9 

62.23531 

3/P005 

144.6 

50.1 

39.9 

57.37138 

4/P013 

144.7 

50.1 

39.9 

51.97498 

5/P009 

144.7 

50.2 

39.8 

45.93798 

6/P003 

144.8 

50.2 

39.8 

39.38992 

7/P008 

144.8 

50.2 

39.8 

32.51828 

8/P010 

144.9 

50.2 

39.8 

25.67801 

9/P012 

145 

50.2 

39.8 

19.3063 

10/P014 

145 

50.3 

39.7 

13.3248 

1 1/P00 1 

145.1 

50.3 

39.7 

8.153563 

12/P015 

145.2 

50.3 

39.7 

3.835664 

13/P011 

145.2 

50.3 

39.7 

-1.84744 

14/P002 

145.3 

50.3 

39.7 

-4.48537 

15/P007 

145.3 

50.4 

39.6 

-7.53443 

Table  6.  Duck,  NC,  Worldview-2  Imagery  data  (March  24,  2010) 
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Name 

Acquisition  Time 
(GMT) 

Mean  Satellite 
Azimuth  (deg) 

Mean  Satellite 
Elevation  (deg) 

1/P002 

19:03:32 

277.8 

61.2 

2/P007 

19:03:42 

267.1 

60.1 

3/P009 

19:03:53 

257.6 

58.2 

4/P006 

19:04:04 

249.5 

55.8 

5/P008 

19:04:15 

242.9 

53.1 

6/P004 

19:04:26 

237.3 

50.1 

7/P010 

19:04:36 

232.9 

47.2 

8/P005 

19:04:47 

229.2 

44.4 

9/P003 

19:04:57 

226.2 

41.6 

10/P001 

19:05:08 

223.6 

39.1 

Name 

Solar 

Azimuth 

(deg) 

Solar 

Elevation 

(deg) 

Solar  Zenith 
Angle  (deg) 

Separation 
Angle  (deg) 

1/P002 

155.9 

56.2 

33.8 

-83.3647 

2/P007 

156 

56.2 

33.8 

-77.9109 

3/P009 

156.1 

56.2 

33.8 

-72.7192 

4/P006 

156.2 

56.3 

33.7 

-68.0675 

5/P008 

156.2 

56.3 

33.7 

-64.2788 

6/P004 

156.3 

56.3 

33.7 

-61.0071 

7/P010 

156.4 

56.3 

33.7 

-58.464 

8/P005 

156.5 

56.3 

33.7 

-56.3613 

9/P003 

156.5 

56.3 

33.7 

-54.8343 

10/P001 

156.6 

56.3 

33.7 

-53.4706 

Table  7.  Pendleton,  CA,  worldview-2  Imagery  data  (March  24,  2010) 
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Table  8.  Satellite  and  Solar  Azimuthal  Angles,  Elevation  Angles,  and  Overall 
Separation  Angles,  Left:  Pendleton,  CA,  Right:  Duck,  NC 
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B.  IMAGE  PREPARATION 

DigitalGlobe  processes  imagery  as  requested  by  the  customer.  Three  levels  are 
available  and  specific  options  may  also  be  available  within  the  level.  These  levels 
include: 

•  Level  1:  Basic  IB  or  Basic  Stereo  Pairs 

•  Level  2:  Standard  (2 A)  or  Ortho-Ready  Standard  (OR2A) 

•  Level  3:  Ortho-rectified  (DigitalGlobe  or  Custom) 

Level  1  processing  is  done  on  all  Level  2  processed  images,  to  include 
radiometric  and  sensor  correction.  Radiometric  correction  is  necessary  to  reduce  the 
effects  of  noise  in  the  imagery.  Noise  can  originate  from  the  atmosphere,  the  sun-sensor 
geometry,  or  the  sensor  itself.  Radiometric  correction  provides  two  key  benefits  to  this 
research.  The  first  allows  mosaicing  of  images.  The  second  advantage  is  image 
consistency,  which  accommodates  work  involving  temporal  data.  The  additional  level  of 
processing  involved  in  Standard  2A  and  OR2A  processing  provides  unifonn  GSD 
throughout  the  image  and  projects  the  image  on  a  plane  using  the  customer’s  selected 
map  projection  and  datum.  Map  projection  options  include  Geographic 
(latitude/longitude),  Universal  Transverse  Mercator  (UTM),  and  a  customer  defined  State 
Plane.  Standard  2A  applies  a  coarse  digital  elevation  model  (DEM),  while  OR2A  is 
projected  to  an  average  elevation,  but  has  no  topographic  relief  applied.  OR2A  stands 
ready  for  custom  ortho-rectification,  providing  flexibility  to  the  customer.  Increased 
processing  comes  at  an  increase  price  to  the  consumer.  With  Level  3  processing,  the 
image  is  ortho-rectified.  Ortho-rectification  applies  a  fine  DEM  using  map  projection 
and  the  customer’s  requested  datum.  The  DEM  can  be  provided  by  DigitalGlobe  or  the 
customer.  Mosaicing  is  also  available  at  the  highest  level  of  processing. 

The  imagery  used  in  this  thesis  underwent  Ortho-Ready  Standard,  a  level  2 
process.  Each  image  over  Duck,  NC,  and  Pendleton,  CA,  was  received  in  three  separate 
tagged  image  fonnat  (TIF)  files.  TIF  files  are  a  standard  imagery  format  as  they  allow 
lossless  imagery  storage  and  can  be  edited  and  resaved  without  degradation  of  the  image. 
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Basic  registration  processes  came  first.  As  the  pass  over  Duck,  NC,  and 
Pendleton,  CA,  resulted  in  15  and  10  images,  respectively,  registration  was  required. 
Registration  is  the  process  of  correcting  a  reference  (warp)  image  to  match  base  image 
geometry,  using  tie  points  between  the  two  images.  While  the  registration  process  is 
tedious,  it  is  necessary.  Errors  arise  due  to  different  satellite  viewpoints  that  cause 
differences  in  acquisition  and  cause  the  images  to  be  misaligned.  In  this  analysis,  the 
chronological-center  image  was  selected  as  the  base  image  to  minimize  the  differences 
between  the  base  image  and  the  warp  images.  While  often  called  ground  control  points 
(GCPs),  the  more  accurate  term  is  tie  points.  Tie  points  refer  to  points  selected  within  the 
base  and  warp  image  and  indicate  a  single  feature  identified  in  both  images.  GCPs  are 
used  to  relate  a  point  in  a  remote  sensing  image  to  a  geographic  point  on  the  earth. 
Registration  allows  the  Geographic  projections  done  on  each  image  to  be  aligned  with  the 
Geographic  projection  of  the  reference  image.  The  ENVI  program  refers  to  tie  points  and 
GCPs  synonymously.  To  achieve  a  well-registered  image,  approximately  70  points  were 
referenced  to  the  base  image.  An  illustration  of  this  process  is  captured  in  Figure  30. 


Figure  30.  The  reference  image  (left)  and  warp  image  (right)  require  the  selection 
of  a  large  number  of  tie  points  for  accurate  registration. 
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When  using  off-nadir  imagery,  the  Geographic  projection  done  by  Level  2 
processing  does  not  result  in  an  accurate  projection  of  latitude  and  longitude,  regardless 
of  the  image  selected  as  the  registration  reference.  If  accurate  latitude  and  longitudes  are 
required,  they  are  accomplished  with  ortho-rectification.  Ortho-rectification  removes  the 
effects  of  terrain  relief  displacement,  lens  distortion,  and  camera  tilt.  Ortho-rectification 
creates  uniform  scale  and  true  geometry  within  imagery  and  allows  measurement  of  true 
distances,  as  an  ortho-rectified  image  accurately  represents  the  Earth's  surface.  This 
correction  is  necessary  when  using  off-nadir  imagery  and  in  topographically  diverse 
landscape.  As  a  satellite  overflies  a  ground  point,  it  images  at  differing  angles  and 
subsequently  returns  different  views  of  the  same  object.  This  occurs  because  as  the  view 
angle  changes,  the  distance  from  the  object  also  changes.  The  same  effect  occurs  when 
the  surface  viewed  is  not  flat,  and  the  effect  becomes  more  pronounced  as  the  diversity  of 
the  terrain  increases.  To  correct  for  terrain  changes,  an  accurate  DEM  is  required.  Figure 
31.  illustrates  how  errors  are  introduced  in  satellite  imagery  due  to  off-nadir  imaging. 
Figure  32.  illustrates  how  terrain  elevation  affects  images. 
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Figure  31. 


(a)  A  flat  surface  viewed  from  directly  overhead  (Nadir)  (b)  The  same 
surface  viewed  with  30°-Off-nadir  Roll  (c)  Surface  with  30  Off-Nadir 
Roll  and  60°  Off-Nadir  Pitch  (d)  Additional  30°  Yaw  Introduced 
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Figure  32.  Displacement  effect  on  imagery  due  to  terrain 

As  the  scope  of  this  thesis  did  not  require  an  accurate  map  of  the  location,  ortho¬ 
rectification  was  not  required.  Registration  of  the  images  was  the  furthest  step  taken  and 
was  done  by  the  staff  of  the  Naval  Postgraduate  School. 

Following  image  registration,  the  images  were  chipped  to  display  smaller  regions. 
Analysis  was  done  on  the  full  image  and  the  smaller  chips. 

C.  REGIONS  OF  INTEREST 

Regions  of  Interest  allow  the  selection  of  an  area  of  pixels  and  allow  them  to  be 
grouped  for  processing  with  ENVI.  ROIs  are  typically  a  group  of  pixels  with  similar 
characteristics.  ROIs  were  selected  for  Duck  and  Pendleton.  To  begin,  a  general  area  of 
terrain  is  selected.  In  Duck,  an  initial  1 1  regions  were  created.  This  eventually  became 
12  as  the  process  revealed  two  types  of  sand  within  one  ROI.  In  Pendleton,  10  regions 
were  selected.  The  ROIs  were  titled  as  follows: 
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Duck,  NC 

Forest 

Pond 

Roof  Top  (x3) 
Baseball  Outfield 
Soil  (Dry) 

Soil  (Wet) 

Ocean 
Golf  Course 
Parking  Lot  (x2) 


Pendleton,  CA 

Field  (Veg) 
Harbor 
Tarmac 
Sand  (x3) 

Soil  (Wet) 

Soil  (Dry) 
Ocean 

Marsh  (Veg) 


Figure  33.  The  initial  1 1  ROIs  at  Duck. 

As  the  ROIs  are  initially  selected  by  a  general  area,  there  are  inherently  some 
dissimilar  pixels  within  the  region.  The  n-D  visualizer  allows  exclusion  of  these  pixels  to 
maintain  a  homogenous  ROI. 
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Figure  34.  ROIs  displayed  in  n-D  visualizer 


The  n-D  visualizer  plots  all  pixels  within  the  ROIs.  Pixel  color  corresponds  with 
the  corresponding  color  of  the  ROI.  As  depicted  in  Figure  34.  ,  the  axes  are  digital 
number  (DN)  in  selected  bands,  with  the  bottom  left  corner  representing  (0,0)  and 
increasing  as  you  move  up  and  right.  In  this  figure,  the  x-axis  is  Band  5  (Red),  and  the  y- 
axis  is  Band  6  (Red-edge).  Note  the  relative  brightness  of  the  vegetation  ROIs  (green 
pixels)  in  the  Red-edge  band  and  the  relative  darkness  in  the  Red  band.  The  tarmac  (pink 
pixels)  is  bright  in  both  bands. 

The  outlier  pixels  are  represented  by  their  separation  from  the  major  grouping  of 
pixels  in  each  ROI.  To  remove  these  from  further  ENVI  processing,  the  user  manually 
highlights  a  group  of  outliers  and  opts  to  remove  them  from  the  group.  Often,  the  n-D 
visualizer  illustrates  clear  differences  between  two  or  more  materials  within  the  same 
group,  which  are  not  easily  discerned  by  the  naked  eye.  ENVI  provides  the  option  of 
splitting  the  ROI  into  two  or  more  new  ROIs,  as  illustrated  with  Figure  35.  When  one 
goes  to  edit  the  Sand  ROI,  Figure  35.  shows  a  pixel  clump.  In  Band  1  and  Band  7,  it  is 
not  very  clear  that  there  are  two  discernible  clumps.  In  Band  5  and  Band  7,  it  becomes 
clearer.  Rather  than  removing  one  of  the  pixel  clumps,  a  second  ROI  is  created.  The 
ROI,  Sand,  is  split  into  Sand  (Dry)  and  Sand  (Wet).  Once  this  process  is  complete,  the 
ROIs  are  updated  accordingly.  The  whittled-down  ROIs  are  displayed  in  Figure  37. 
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Figure  35.  The  n-D  visualizer  display.  Left:  Sand  ROI  in  (Band  1  &  7),  Center: 

Two  identified  clumps  (Bands  1  &  7),  Right:  Two  clumps  in  Band  5  &  7, 
which  introduces  more  separate  in  the  clumps 


Figure  36.  Left:  Sandy  area,  Center:  Initial  Sand  ROI  which  encompasses  wet  and 
dry  sand,  Right:  The  result  of  n-D  visualizer.  The  original  Sand  ROI  is 
divided  into  2  different  ROIs,  wet  and  dry. 


These  pixels  typically  represent  differing  materials  within  the  terrain  such  as  a 
road  cutting  through  a  field  or  a  car  in  a  parking  lot.  By  removing  these  from  processing, 
a  homogenous  ROI  can  be  obtained  for  analysis. 

While  at  DigitalGlobe,  the  product  underwent  processing  that  made  use  of  the 
gain  and  offset  numbers  provided  in  the  Appendix.  These  numbers  were  retrieved  from 
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the  .imd  files  provided  along  with  the  imagery.  The  end  result  delivered  to  the  customer 
transfonned  the  raw  image  digital  number  values  into  radiometrically  corrected  image 
pixels,  qpixei,band-  The  calibration  factor  applied  can  also  be  found  in  the  Appendix. 

Statistics  were  run  on  the  ROIs.  Worldview-2  DNs  were  converted  to  absolute 
radiance  utilizing  the  following  equations: 

DN  =  q pad,  band  [counts]  (18) 

Lpixei,  band  =  abscaIfactorba„d  *  q pixel,  band  [Wm  2sr  1  ]  (19) 


LXptxel,  band  — 


Lpixei,  band  r  _2  .  ] 


AXband 


[Wm  sr  pm] 


(20) 


where  qpixei,band  are  radiometrically  corrected  image  pixels,  Lpixci_hand  is  top-of-atmosphere 
band  integrated  radiance  image  pixels,  and  Lx,pixei,band  is  top-of-atmosphere  band-averaged 
spectral  radiance  image  pixels.  Following  this,  statistical  analyses  were  then  done  on  the 
homogeneous  ROIs  utilizing  ENVI,  and  the  radiance  values  were  analyzed. 


Figure  37.  (Left)  Several  of  Duck,  NC,  original  ROIs.  (Right)  n-D  visualizer  has 
removed  the  pixels  that  displayed  differing  DNs  from  the  majority  over 
several  considered  bands,  representing  cars  in  a  parking  lot  (Red  ROI),  an 
island  in  a  pond  (Blue  ROI)  and  structures/objects  on  a  roof  (Pink  ROIs) 
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D.  IMAGE  RATIOING 


Image  ratioing  offers  the  benefit  of  removing  or  minimizing  many  of  the  errors 
previously  mentioned.  Image  ratioing  is  the  process  of  dividing  an  image  by  another 
image,  pixel  by  pixel.  The  process  was  done  using  the  ENVI-compatible  IDL  workbench 
and  was  performed  across  all  inputted  wavelengths. 


Figure  38.  Image  Ratioing  Interactive  Data  Language  Code 

In  image  ratioing,  the  brighter  image  is  divided  into  the  darker  image  to  ensure 
the  resulting  image  is  scaled  so  that  all  pixels  hold  a  value  between  0  and  1.  This  was 
done  for  both  sets  of  imagery.  In  the  Pendleton  imagery,  image  P009  was  divided  by 
P010,  and  in  the  Duck  imagery  P006  was  divided  by  P011.  In  both  cases,  the  brightest 
image,  placed  in  the  denominator,  corresponded  to  the  up-sun,  or  on-axis,  image  where 
the  satellite  elevation  most  closely  matched  the  solar  elevation.  Coincidentally,  the 
brightest  image  was  the  latter  image  taken.  The  resulting  ratio  represents  congruity 
between  the  initial  two  images.  Zero  represents  dissimilar  pixels,  such  as  a  saturated 
pixel  contrasted  with  a  dark  pixel,  while  a  value  of  1  represents  identical  pixels.  Values 
lying  between  represent  some  disparity  in  pixel  value.  This  concept  is  illustrated  in 
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Figure  39.  Note  the  boat  traveling  south  down  the  coast.  It  appears  bright  in  contrast 
with  the  surrounding  ocean-  a  high  DN  value  in  contrast  with  a  low  DN.  As  mentioned 
previously,  the  brighter  image  happened  to  be  the  second  image  taken  and  was  used  in 
the  denominator.  At  the  boat’s  position  in  the  first  image,  its  high  value  is  placed  in  the 
numerator.  At  that  same  location  in  the  second  image,  there  is  only  dark-pixelated  ocean, 
which  is  placed  in  the  denominator.  This  results  in  a  value  near  1.  The  values  are 
reversed  at  the  boat’s  second  position,  resulting  in  a  value  near  0. 


Figure  39.  Image  Ratio  of  Duck,  NC,  Image  P006/P01 1  (Band  2:  450-510  nm) 


The  grayscale  image  displays  information  for  only  one  band.  Figure  40.  displays 
information  for  three  bands  and  the  use  of  color  in  the  image  allows  for  the  presentation 
of  a  greater  depth  of  information.  As  in  the  grayscale  image,  white  and  black  represent 
total  similarity  and  total  dissimilarity,  respectively;  however,  in  the  color  image  the  range 
of  gray  now  displays  equal  change  between  the  three  represented  bands.  Gray  is  useful 
for  displaying  the  level  of  view-angle  dependence  of  reflectance.  Color  is  used  to 
illustrate  unequal  change  between  the  bands  and  illustrates  both  view-angle  dependence 
and  wavelength  dependence  of  reflectance.  The  color  image  of  Duck,  NC,  is  displayed 
in  Figure  40. 
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Figure  40.  Image  Ratio  of  Duck,  NC,  Image  P006/P01 1  (Bands  2:  450-510,  Band 

5:  630-690,  and  Band  8:  860-1040  mn) 


Image  ratioing,  in  addition  to  allowing  the  presentation  of  information  in  a 
concentrated  form,  offers  the  benefit  of  a  reduction  in  severity  of  numerous  undesirable 
effects.  The  images  were  taken  in  close  sequential  order,  with  the  solar  zenith  angle 
nearly  identical;  therefore,  there  is  similar  shadowing  between  the  images.  Ratioing 
removes  a  large  amount  of  this  shadowing,  and  presents  the  shadowed  regions  in  a 
manner  that  is  more  consistent  with  the  sunlit  regions.  Additionally,  the  satellite  views 
the  target  through  the  Earth’s  atmosphere.  The  atmosphere  has  major  effects  on  imagery, 
and  not  all  affect  different  wavelengths  equally.  Scattering  is  one  such  affect,  which 
causes  shorter  wavelengths  to  be  scattered  more  than  longer  wavelengths.  For  example, 
the  majority  of  atmospheric  effects  are  cancelled  out  by  ratioing.  Some  effects  remain 
due  to  the  satellite’s  motion,  which  causes  the  view  angle,  and  subsequently  the  path 
length,  to  change.  Furthermore,  in  our  analysis,  accurate  in  situ  solar  irradiation 
measurements  were  not  taken.  These  measurements  are  part  of  the  BRDF  calculation, 
and  their  removal  through  ratioing  allows  the  comparison  of  radiance  to  more  accurately 
scale  to  the  BRDF.  Analysis  on  these  images  is  done  in  the  next  chapter. 
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V.  OBSERVATIONS  AND  ANALYSIS 


A.  CALCULATIONS 

Calculations  were  done  using  the  information  regarding  the  satellite’s  and  the 
sun’s  positions  when  the  imagery  was  acquired.  The  separation  angle  between  the 
incident  angle  (the  sun’s  position)  and  the  reflected  angle  (the  satellite’s  position)  were 
calculated  to  allow  comparison  between  images.  This  information  is  available  in  the 
Appendix. 


Figure  41 .  Illustration  of  Solar  Zenith  Angle,  Viewing  Zenith  Angle,  and 

Calculated  Separation  Angle  [36] 


B.  EFFECTS  OF  LIMITED  IMAGERY  PROCESSING 

The  raw  imagery  used  in  this  research  is  subject  to  many  effects  that  must  be 
accounted  for  in  order  to  derive  quantitative  data.  Limited  processing  was  performed  on 
the  imagery,  and  many  of  these  effects  remain.  Processing  limitations  must  be  accurately 
annotated  to  ensure  inaccurate  analysis  is  not  obtained  or  assumed.  Initial  processing  was 
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done  by  DigitalGlobe  prior  to  its  dissemination  to  the  customer.  Upon  receipt  of  the 
imagery  at  the  Naval  Postgraduate  School,  all  image  processing  was  done  utilizing  the 
ENVI  software.  Processing  tools  used  consisted  of  the  Worldview  2  Calibrated  Utility 
preprocessing  tool,  Georeferenced  mosaicing,  and  Image-to-Image  registration.  No 
further  processing  tools  were  applied  to  the  imagery  prior  to  follow-on  analysis. 
Therefore,  it  can  be  assumed  that  all  effects  not  corrected  for  by  these  tools,  remain. 
Specifically,  effects  imparted  due  to  atmospheric  conditions  are  not  removed. 
Additionally,  orthorectification  was  not  performed  on  the  imagery,  and  the  spatial 
processing  done,  registration,  does  not  accurately  allow  for  accurate  DEM  integration. 

C.  MISR  DATA 

A  comparable  data  set  was  obtained  for  Pendleton  from  MISR.  MISR  data  were 
requested  and  received  from  the  EOS  data  gateway  at  NASA  Langley’s  Atmospheric 
Science  and  Data  Center  (ASDC).  As  MISR  images  a  geographic  area  approximately 
once  every  16  days,  the  imagery  requested  was  taken  on  March  26,  2010.  This 
corresponded  to  Terra  orbit  54,630  and  path  41,  as  depicted  in  Figure  42.  This  data  was 
selected  because  it  was  the  timeliest  for  comparison  with  the  imagery  taken  by 
Worldview-2  on  March  24,  2010. 

The  product  received  was  Level  1B1  Radiance  Data  in  Global  mode  for  each  of 
MISR’s  9  cameras.  The  level  of  processing  included  DNs  radiometrically-scaled  to 
radiance  with  no  geometric  resampling  or  registration.  Following  receipt,  the  product 
was  processed  using  ENVI  and  chipped  to  the  desired  region.  The  image  was  relatively 
cloud  free  in  the  region  of  interest. 
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Figure  42. 


Figure  43. 


Pendleton  Overflight  by  Terra/MISR  (Path  41)  on  the  WRS-2  Map 
(Worldwide  Reference  System) 


MISR  Radiance  Image:  ROI  for  comparison  with  Worldview-2 
(Camera:  AN,  Red  Band) 
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Comparison  between  the  Worldview-2  and  MISR  were  done  utilizing  a  ROI 
generated  in  the  coastal  hills  in  the  vicinity  of  Pendleton.  While  the  spatial  size  of  the 
image  was  similar,  approximately  125  km  ,  the  spatial  resolution  of  each  image  resulted 
in  a  marked  difference  in  pixels  for  analysis.  In  global  mode,  spatial  resolution  for  MISR 
is  approximately  1  km  for  all  channels  except  the  Red  channel,  which  is  275  m,  while 
Worldview-2’s  spatial  resolution  is  2.4  m  (off-nadir).  Analysis  was  done  on  MISR’s 
scaled  red  channel  and  Band  5  of  Worldview-2.  The  MISR  red  channel  was  selected  due 
to  its  high  spatial  resolution  in  contrast  with  the  other  three  bands. 

While  the  MISR  product  returned  scaled  radiance  values,  the  Worldview-2  data 
underwent  conversion  from  DN  to  radiance  utilizing  the  process  describe  in  the  Region 
of  Interest  Section  of  this  chapter.  The  trends  are  presented  for  comparison. 

Figure  46.  shows  a  decrease  in  radiance  as  the  zenith  angle  decreases.  This  is  in 
contrast  with  measurements  made  by  systems  such  as  Chen-Leblanc’s 
gonioreflectometer,  PARABOLA,  and  POLDER  (Figure  13.  ,  Figure  15.  ,  Figure  21.  , 
respectively).  While  most  plots  show  a  comparison  of  the  incident  and  viewing  angles, 
an  estimate  of  the  sun’s  position  and  the  corresponding  separation  angle  were  made.  This 
does  not  appear  to  provide  an  explanation  for  the  difference;  however,  the  data  does  not 
appear  to  be  erroneous  as  the  right  half  of  Figure  47.  displays  a  similar  trend  to  the  MISR 
data  in  Figure  46.  The  cause  for  the  difference  in  the  region  is  unknown;  however,  the 
region  is  large  and  diverse,  and  it  can  be  assumed  that  the  materials  present  in  the  region 
are  varied.  As  illustrated  in  the  PROBA/CHRIS  data  of  Figure  23.  ,  different  materials 
display  different  angle-dependency  trends,  and  this  is  one  possible  explanation  for  the 
unexpected  radiance-vs. -angle  curves. 
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Figure  44.  MISR  Radiance  Image:  Zoom  of  ROI  for  comparison  with  Worldview -2 

(Camera:  AN,  Red  Band) 


Figure  45.  Worldview-2  Image:  ROI  for  Comparison  with  MISR  (Image  P009, 

Band  5:  630-690  mn) 
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MISR:  26  MAR  2010 


Zenith  Angle  (deg) 


Figure  46.  MISR  Radiance  Plot  in  Red  Channel  (672  ran)  by  Camera  and  Zenith 

Angle:  Pendleton,  CA 


Worldview-2:  24  MAR  2010 


Figure  47.  Worldview-2  Radiance  Plot  in  Band  5  (630-690  nm)  by  Time  and 

Zenith  Angle:  Pendleton,  CA 
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Occurences 


Figure  48.  MISR  Brightness  Variance  in  Red  Channel  (672  rnn)  with  respect  to 

Elevation  Angle  (Camera):  Pendleton,  CA 


endleton,  CA  -  Band 


Figure  49.  Worldview-2  Brightness  Variance  in  Band  5  (630-690  nm)  with  respect 

to  Elevation  Angle  (Time):  Pendleton,  CA 
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It  is  difficult  to  accurately  compare  individual  ROIs  given  the  enormous  spatial 
scale  difference  between  the  systems.  We  approach  the  data  globally,  therefore,  by  using 
scene  histograms,  as  shown  in  Figure  48.  and  Figure  49.  Figure  48.  shows  how  the  672 
nm  MISR  data  vary  for  the  image  chip  shown  in  Figure  43.  as  the  view  angle  varies.  The 
nadir  view  has  the  dark  peak,  the  most  oblique  views  are  the  brightest.  The  peak  is 
largely  the  dark  water  which  dominates  the  scene.  The  ROI  was  evaluated  in  the  red 
channel  of  each  of  the  nine  images,  and  the  radiance  values  were  plotted.  Radiance 
shows  a  decrease  with  increasing  elevation  angle.  This  is  also  shown  in  the  whole  image 
histogram  plotted  in  Figure  48.  In  the  latter  figure,  the  radiance  curve  shows  a  leftward 
shift  until  reaching  nadir,  after  which  it  shifts  back  to  the  right.  The  images  were 
acquired  around  noon  over  Pendleton,  CA.  Figure  49.  shows  the  same  sort  of  shift  in  the 
‘water’  peak,  but  also  a  contraction  in  the  dynamic  range,  which  is  so  far  unexplained. 

D.  SHIFT  IN  DIGITAL  NUMBER  OBSERVED  OVER  SEQUENCE  OF 

IMAGES 

Analysis  was  done  on  the  imagery  of  Duck,  NC,  and  Pendleton,  CA.  The  large 
spike  in  the  data  in  the  low  DNs  is  the  water  in  the  scene.  The  ocean  is  a  poor  reflector  of 
energy,  especially  in  the  NIR  bands.  Conversely,  plants  are  excellent  NIR  reflectors  and 
this  is  captured  in  the  spike  at  the  DN  1400  count.  The  tail  of  the  curve  (>1400  counts) 
represents  the  extremely  reflective  nature  of  man-made  materials  within  the  image. 

A  brightness  shift  across  all  bands  was  observed.  This  observed  shift,  represented 
in  the  figure,  qualitatively  represents  the  BRDF  of  the  scene.  As  the  images  progress 
chronologically,  represented  by  an  inverse  rainbow  scheme,  there  is  an  observed  shift  to 
the  right  of  DN  through  the  first  13  images,  followed  by  a  shift  back  to  the  left  for  the  last 
2  images.  Graphs  of  all  eight  bands  for  both  Duck  and  Pendleton  can  be  found  in  the 
Appendix. 
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Duck  NC 


Bond 


Figure  50.  Worldview-2  Brightness  Variance  in  Band  8  (860-1040  inn)  with 

respect  to  Elevation  Angle  (Time):  Duck,  NC 


E.  IMAGE  RATIO  ANALYSIS 
1.  Images 

Figure  51.  is  a  full  image  ratio  of  Duck,  NC,  in  the  Worldview-2  bands  2,  5,  and 
8.  Figure  52.  is  a  full  image  ratio  of  the  Pendleton  scene  in  the  same  bands.  While  the 
images  acquired  over  Duck  appear  relatively  cloud-free,  there  is  a  large  amount  of  cloud 
cover  within  the  Pendleton  images,  which  affects  the  BRDF  values. 
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Figure  51.  Duck  Image  Ratio,  Image  P006/P01 1  (Bands  2:  450-510,  Band  5:  630 

690,  and  Band  8:  860-1040  mn) 


Pendleton  Image  Ratio,  Image  P009/P010  (Bands  2:  450-510,  Band  5 
630-690,  and  Band  8:  860-1040  mn) 
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Figure  52. 


The  full  Pendleton  scene  was  then  chipped  down  to  display  interesting  regions,  as 
is  the  case  in  Figure  53.  and  Figure  54.  These  figures,  as  well  as  Figure  55.  ,  display  how 
BRDF  varies  over  view-angle  and  wavelength  despite  similar  terrain.  Three  fields  of 
differing  make  up  are  viewed  in  the  figures.  In  the  first  image,  a  gray-scale  color  ramp 
shows  the  range  of  ratio  values  within  the  image  in  Band  3.  Darker  regions  indicate  a 
greater  difference  between  the  compared  images,  and  therefore  a  greater  dependence  on 
view  angle.  In  the  second  image,  the  same  concept  is  captured  through  the  representation 
of  differing  colors;  however,  the  image  also  illustrates  the  BRDF  dependence  on 
wavelength. 


Figure  53.  Pendleton  Image  Ratio:  ROIs  annotated  (Band  3:  510-580  nm) 


Further  analysis  of  this  was  done  on  selected  ROIs.  A  better  understanding  of  the 
colors  observed  in  Figure  54.  can  be  achieved  through  observing  the  ratio  changes 
annotated  in  Figure  55.  As  none  of  the  regions  have  a  flat  line  on  the  relative  ratio  chart, 
each  region  has  an  element  of  color.  The  wet  and  dry  soils  have  similar  wavelength- 
dependent  BRDF  shifts  across  all  wavelengths,  represented  by  the  two  curves 
approximately  paralleling  each  other.  If  the  ratios  moved  in  lock  step,  the  areas  would 
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differ  only  in  the  brightness  of  color.  They  differ  slightly  in  color,  however,  because  the 
curves  differ  more  than  by  solely  a  y-axis  displacement.  As  it  can  be  observed  in  Figure 
55.  ,  the  two  regions  are  more  similar  than  the  third  region,  the  Field  (Veg).  The  relative 
ratio  of  the  vegetated  field  ticks  up  closer  to  1  in  the  longer  wavelengths.  This  allows  one 
to  conclude  that  vegetation  has  a  reduction  in  angular  dependence  in  these  wavelengths. 
The  marked  shift  between  the  bands  results  in  a  brightly  colored  region.  All  3  ROIs 
move  from  the  bottom  left  to  the  top  right.  This  indicates  that  Band  8  appears  to  be  the 
least  affected  by  changes  in  the  view  angle. 


Figure  54.  Chip  of  Pendleton  Image  Ratio:  Field  (Veg)  and  Soil  ROIs  annotated 
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Figure  55.  Relative  Ratio  of  Pendleton  Image  Ratio:  Field  (Veg)  and  Soil  ROIs 
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2. 


Terrain  Effects  in  Spatial  Profiles 


Figure  56.  Google  maps  true-color  image  of  Pendleton  Spatial  Profiles’  Regions 


Figure  57.  Pendleton  Image  Ratio  Spatial  Profiles  (Band  3:  510-580  nm) 
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Figure  56.  and  Figure  57.  are  images  of  an  area  in  the  Pendleton  scene.  Figure  56. 
is  a  true-color  image  acquired  through  Google  maps,  while  Figure  57.  is  a  chip  from  the 
image  ratio  of  P009  and  P010  in  Band  3.  Spatial  profiles  were  created  within  these 
images.  They  sample  the  value  of  pixels  along  the  line.  They  were  drawn  with  the  intent 
of  spanning  a  region  of  changing  elevation  and  terrain.  Profile  1  spans  a  valley,  and  this 
valley  can  be  seen  in  the  approximate  locations  of  50-100  in  Figure  58. 
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Figure  58.  Pendleton  Image  Ratio  Spatial  Profile  #1  (Band  3:  510-580  nm) 


Profile  2  is  across  a  mountain  range,  with  the  left  side,  0-40,  corresponding  to  a 
westerly  facing  slope,  and  the  right  side  facing  east,  displayed  in  Figure  59.  These 
profiles  represent  how  BRDF  changes  while  the  ground  angle  changes,  therefore,  the 
exitance  angle  is  changing,  but  the  viewing  angle  remains  fixed. 


Spatial  Profile  1:  Band  3 
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Figure  59.  Pendleton  Image  Ratio  Spatial  Profile  #2  (Band  3:  510-580  nm) 


F.  BRDF  SHIFT  BETWEEN  REGIONS  OF  INTEREST  (DUCK,  NC) 

The  sun  and  the  satellite’s  position  relative  to  the  sun  are  plotted  in  Figure  60. 
The  series  of  separation  angles  are  depicted  with  an  arc  with  plus  marks  at  each  specific 
angle.  These  plus  marks  are  colored  to  reflect  the  relative  BRDF  value  and  the  color 
ramp  provides  the  key. 

Wet  and  dry  soil  in  band  7,  the  NIR  band,  is  compared  in  Figure  60.  The 
materials  are  similar,  therefore  the  differences  are  subtle,  but  it  can  be  observed  that  the 
dry  soil  exhibits  a  generally  higher  BRDF  for  all  plotted  angles.  This  is  consistent  with 
general  spectral  measurements  of  soil.  Soil  reflectivity  is  largely  dependent  on  six 
variables:  moisture  content,  organic  matter  content,  particle  size  distribution,  iron  oxide 
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content,  soil  mineralogy,  and  soil  structure,  with  moisture  content  having  the  largest 
impact.  This  suggests  the  difference  between  the  two  is  caused  by  the  additional 
absorption  by  the  water  in  the  wet  soil. 


BRDF 

Scaled  *  100 


BRDF 
Scaled*  100 


Figure  60.  BRDF  ROI  Comparison  wrt  Satellite  Position  for  Duck  ROIs: 

Left:  Soil  (Wet),  Right:  Soil  (Dry)  (Band  7:  770-895  nm) 


G.  BRDF  SHIFT  BETWEEN  BANDS 


Figure  6 1 .  BRDF  Band  Comparison  wrt  Satellite  Position  for  Duck  ROI: 

Soil  (Dry)  Band  7  and  8  (Band  7:  770-895  nm,  Band  8:  860-1040  nm) 
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Another  set  of  BRDFs  are  compared  in  Figure  61.  This  time  the  difference  in 
BRDF  across  different  bands  is  explored.  Dry  soil  is  plotted  in  the  figure  with  band  7 
BRDFs  on  the  left  and  band  8  BRDFs  on  the  right.  The  figure  on  the  right  displays  a 
higher  BRDF  value  across  all  separation  angles.  This  is  also  consistent  with  soil  studies, 
which  show  reflectance  increasing  until  the  SWIR  band  [37]. 

H.  IMAGE  RATIO  DISPLAYING  WAVELENGTH  DEPENDENCE  IN 
REGIONS  OF  INTEREST  BY  BAND 


The  plot  of  image  ratio  relative  ratios  for  Duck,  NC,  and  Pendleton,  CA,  are 
depicted  in  Figure  62.  and  Figure  63.  ,  respectively. 
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Figure  63.  Image  Ratio  done  for  All  Regions  of  Interest  in  all  Worldview -2  Bands 

(Pendleton,  CA) 

Numerous  material  characteristics  can  be  extracted  from  the  charts.  ROI  curves 
that  are  flat  indicate  that  there  is  not  a  large  dependency  on  wavelength  (Duck:  Parking 
Lot,  Pendleton:  Sand  1,  2,  and  3).  Flat  curves  close  to  1  also  indicate  a  low  dependency 
on  angle.  If  the  curve  moves  dramatically  across  wavelengths,  this  reflects  a  wavelength 
dependency  of  the  BRDF  (Duck:  Baseball  Field,  Pendleton:  Field  (Veg)).  Ratio  values 
approaching  a  value  of  1  are  indicative  of  similarity  between  the  two  input  images  and 
low  dependence  on  angle  (Duck:  Roof  Top),  while  low-value  ratios  convey  that  the  input 
images  are  dissimilar  in  the  respective  band  due  to  the  viewing  angle  (Duck:  Forest). 
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I. 


BRDF  RESULTS  BY  SEPARATION  ANGLE  FOR  REGIONS  OF 
INTEREST 


When  displayed  by  separation  angle,  a  clear  pattern  in  the  radiance  and 
subsequently  the  BRDF  is  observed.  Studying  the  charts  in  the  appendix  reveals 
information  about  the  surfaces  within  the  ROIs.  No  surfaces  in  nature  are  a  perfect 
example  of  a  specular  or  diffuse  reflector  and  are  instead  less-idealized  versions  of  the 
two.  A  less-idealized  surface  may  appear  brighter  in  the  specular  and  backscatter 
directions  and  darker  when  viewed  at  grazing  angles.  The  choppiness  of  water  affects  its 
reflective  characteristics. 


Figure  64.  BRDF  Results  by  Separation  Angle  for  Ocean  vs.  Harbor  (Pendleton, 

CA) 


As  observed  in  Figure  64.  ,  the  ocean  is  presumably  rougher  than  the  inland 
harbor  and  therefore  it  exhibits  less  specular  traits.  The  BRDF  of  the  harbor  rises  in  the 
specular  direction  and  again  as  the  angle  decreases  toward  the  backscatter  direction. 
Spectral  trends  can  be  observed,  for  example,  water  reflects  similarly  in  both  the  near-IR 
and  the  red  bands,  seen  in  Figure  64.  Another  spectral  effect  is  seen  in  Figure  65. 
Vegetation  is  highly  reflective  in  the  IR  and  absorbs  strongly  in  the  red,  which  can  be 
clearly  seen  in  the  figure  Band  5  and  Band  8  curves.  Vegetation  is  an  example  of  a  rough 
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surface  and  typically  has  a  strong  return  from  the  backscatter  direction,  called  the  hot 
spot.  This  effect  can  be  seen  in  Figure  65.  It  does  not  show  a  specular  component  in  the 
forward-scatter  direction,  but  rises  as  the  separation  angle  decreases. 


Figure  65.  BRDF  Results  by  Separation  Angle  for  Field  (Veg)  (Pendleton,  CA) 


Figure  66.  displays  the  highly  reflective  nature  of  man-made  materials.  These 
materials  are  reflective  across  all  wavelengths,  and  display  nearly  diffuse  characteristics, 
corresponding  to  a  high  BRDF  in  all  wavelengths  and  across  all  angles. 


Duck:  Roof  Top:  BRDF 


Pencleton:  Tarmac:  BRDF 


Figure  66.  BRDF  Results  by  Separation  Angle  for  man-made  materials 
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Some  materials  reflect  IR  energy  better  than  others.  This  can  be  seen  in  Figure 
67.  Water  is  inherently  a  poor  reflector,  and  this  is  particularly  true  in  the  near-IR  and  IR 
bands.  Vegetation  exhibits  the  opposite  effect  in  the  near-IR  and  IR. 
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Figure  67.  BRDF  Results  by  Separation  Angle  for  Excellent  and  Poor  IR 

Reflectors 
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VI.  SUMMARY 


Bidirectional  Reflectance  Distribution  Function  (BRDF)  trends  were  derived  from 
the  8  spectral  bands  of  Worldview-2  imagery,  taken  over  Duck,  NC,  and  Pendleton,  CA, 
on  March  24,  2010.  Twenty-five  Level  OR2A  images  were  retrieved  from  DigitalGlobe, 
after  which  they  underwent  co-registration  and  processing  utilizing  the  ENVI  software. 

ENVI  tools  provided  the  means  to  analyze  the  imagery.  Histograms  of  the  images 
were  generated.  From  these  charts,  BRDF  shifts  were  visually  apparent  and  could  be 
compared.  Analysis  of  the  data  also  included  the  use  of  ENVI’s  ROI  tool  and  the  n-D 
visualize  tool.  With  these  tools,  numerous  areas  were  highlighted  and  refined  for 
examination.  These  regions  represented  a  variety  of  diverse  terrains,  to  include: 
vegetation,  asphalt,  sand,  and  water.  Image  ratioing  was  completed  on  select  images 
utilizing  the  compatible  IDL  program,  which  allowed  the  removal  of  many  undesired 
imagery  effects.  Spatial  profiles  were  generated  from  the  created  image  ratios.  ENVI’s 
pre-processing  tool  allowed  the  conversion  of  DN  to  radiance  values,  which  scaled  the 
observed  values  to  allow  reasonable  BRDF  values  to  be  generated;  however,  the  values 
contained  within  this  document  are  not  independently  meaningful,  but  instead  the  trends 
they  display  and  the  methods  of  retrieval  examined  are  valuable. 

Separation  angles  between  the  inbound,  solar  zenith  angle  and  the  outbound 
reflectance  angle  were  calculated  for  comparison.  Plotted  against  these  numbers,  BRDF 
shifts  across  wavelengths  and  regions  of  interest  were  examined. 

Data  was  pulled  from  NASA  Langley’s  ASDC  from  March  26,  2010  for  a  rough 
comparison  over  Pendleton,  CA.  Level  1B1  Radiance  products  from  MISR’s  9  multi¬ 
angle  cameras  were  downloaded  and  chipped  for  further  evaluation  with  ENVI.  Within 
the  MISR  data  and  the  Worldview-2  data,  a  common  ROI  was  generated  and  the 
accompanying  statistics  scrutinized. 
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VII.  CONCLUSIONS  AND  FUTURE  WORK 


Worldview-2  is  an  important  tool  in  the  expanding  field  of  global  bidirectional 
reflectance  distribution  function  mapping,  with  applications  in  areas  like  terrain  mapping 
[38]. 

Worldview -2  provides  spatial  resolution  at  a  level  that  had  previously  not  been 
obtainable  by  on-orbit  remote  sensors.  Aerial  photography  allows  for  a  similar  spatial 
resolution,  however,  satellites  provide  high  revisit  times  that  would  be  cost  prohibitive  to 
the  aerial  imagery  community.  Worldview-2  provides  images  at  a  spectral  resolution  that 
is  beneficial  to  the  BRDF  study.  While  Worldview-2’s  multispectral  imager  does  not 
offer  the  large  set  of  spectral  data  that  a  hyperspectral  imager  would,  its  use  is 
advantageous.  On-orbit  hyperspectral  sensors  are  disadvantageous  as  they  increase 
processing  and  increase  the  complexity  and  cost  of  the  satellite’s  operation;  however, 
DigitalGlobe’s  judicious  selection  of  8  bands  provides  access  to  beneficial  data  at  key 
wavelengths.  Based  on  the  results  of  numerous  aerial  hyperspectral  studies,  8  spectral 
bands  at  wavelengths  critical  in  the  study  of  change  detection  were  selected. 

This  thesis  delineates  a  technique  to  further  the  growth  of  the  library  of  BRDF 
measurements,  and  utilizes  tools  capable  of  extracting  these  measurements. 

Future  work  on  this  subject  should  first  focus  on  improving  the  process  of 
extracting  reflectance  data  from  satellite  imagery.  ENVI  provides  an  ortho-rectification 
tool,  which  was  not  utilized  in  this  thesis;  exploring  its  use  could  prove  beneficial. 
Imagery  was  obtained  from  multiple  space-borne  sensors  but  was  not  coordinated  to  be 
obtained  simultaneously.  Advanced  coordination  between  ground  and  space-borne  assets 
could  provide  corroboration  for  data.  Processing  the  data  was  time  intensive  and 
therefore  this  thesis  utilized  only  25  data  points.  To  obtain  a  comprehensive  view  of  the 
reflectance  environment,  it  is  necessary  to  model  the  BRDF  with  significantly  more  data 
points.  With  the  more  expeditious  processing  technique  of  ortho-rectification, 
substantially  more  data  could  be  handled.  With  the  process  improved,  a  real  contribution 
of  data  of  high-spatial  resolution  in  8-bands  could  be  submitted  to  the  BRDF  library. 
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APPENDIX:  CHARTS  AND  GRAPHS 


A.  WORLDVIEW-2  GAIN,  OFFSET  VALUES*,  ABSOLUTE  CALIBRATION 
FACTOR  AND  EFFECTIVE  BANDWIDTH 


Band 

Gain 

Absolute  Calibration 

Factor 

Effective  Bandwidth  (pm) 

1 

1.85415428e-002 

9.295654e-003 

4.73e-002 

2 

9.81 1 18344e-003 

1.260825e-002 

5.43e-002 

3 

1.55439861e-002 

9.713071e-003 

6.30e-002 

4 

1.07401758e-002 

5.101088e-003 

3.74e-002 

5 

1.07401758e-002 

1.103623e-002 

5.74e-002 

6 

1.07401758e-002 

4.539619e-003 

3.93e-002 

7 

1.07401758e-002 

1.224380e-002 

9.89e-002 

8 

1.07401758e-002 

9.042234e-003 

9.96e-002 

*  Offset  for  all  bands  equaled  zero 


B.  SPECTRAL  IRRADIANCE 


Band 

Spectral  Irradiance 
(W/(m2*sr*nm)) 

1 

1.88489 

2 

2.05999 

3 

2.09716 

4 

2.00309 

5 

1.8758 

6 

1.68533 

7 

1.36722 

8 

1.16124 

^Generated  by  www.spectralcalc.com  (T  =  5250  C,  s=l) 
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c. 


MISR  PRODUCTS  UTILIZED 


Type 

Identifier 

Start  Time  (GMT) 

End  Time  (GMT) 

MI1B1 

MISRAM  1RPGMP04  l_0054630_DF_F03_0024.hdf 

18:25:12.000000 

19:15:09.160000 

MI1B1 

MISRAM  1RPGMP04  l_0054630_CF_F03_0024.hdf 

18:26:13.440000 

19:16:10.600000 

MI1B1 

MISRAM  1RPGMP04  l_0054630_BF_F03_0024.hdf 

18:27:07.712000 

19:17:04.872000 

MI1B1 

MISRAM  1RPGMP04  l_0054630_AF_F03_0024.hdf 

18:27:54.816000 

19:17:51.976000 

MI1B1 

MISR_AM  1_RP_GM_P04  l_0054630_AN_F03_0024.hdf 

18:28:41.920000 

19:18:39.080000 

MI1B1 

MISRAM  1RPGMP04  l_0054630_AA_F03_0024.hdf 

18:29:29.024000 

19:19:26.184000 

MI1B1 

MISRAM  1RPGMP04  l_0054630_BA_F03_0024.hdf 

18:30:16.128000 

19:20:13.288000 

MI1B1 

MISRAM  1RPGMP04  l_0054630_CA_F03_0024.hdf 

18:31:10.400000 

19:21:07.560000 

MI1B1 

MISRAM  1RPGMP04  l_0054630_DA_F03_0024.hdf 

18:32:11.840000 

19:22:09.000000 

*Data  acquired  on  March  26,  2010.  Provided  by  EOS  data  gateway  at  NASA  Langley’s  Atmospheric  Science  and  Data  Center. 


Key 

MI  =  MISR 

1B1  =  Level  1B1  Product 

MISR  =  Instrument  Name 

AMI  =  Satellite  Name 

(G)RP  =  (Georectified)  Radiance  Product 

GM  =  Global  Mode 

Pxxx  =  Path  Number 

Oxxxxxx  =  Orbit  Number 

XX  =  Camera 

Lnn  =  Lormat  version  of  the  product 
nnnn  =  Product  Version  number 
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D.  HISTOGRAM  SHIFT  WITH  VIEW  ANGLE  BY  BAND  (DUCK,  NC) 
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Figure  68.  Histogram  Shift  with  View  Angle:  Duck,  NC  (Band  1) 
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Figure  69.  Histogram  Shift  with  View  Angle:  Duck,  NC  (Band  2) 
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Figure  70.  Histogram  Shift  with  View  Angle:  Duck,  NC  (Band  3) 
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Figure  7 1 .  Histogram  Shift  with  View  Angle:  Duck,  NC  (Band  4) 
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Figure  72.  Histogram  Shift  with  View  Angle:  Duck,  NC  (Band  5) 
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Figure  73.  Histogram  Shift  with  View  Angle:  Duck,  NC  (Band  6) 
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Figure  74.  Histogram  Shift  with  View  Angle:  Duck,  NC  (Band  7) 
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Figure  75.  Histogram  Shift  with  View  Angle:  Duck,  NC  (Band  8) 
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E.  HISTOGRAM  SHIFT  WITH  VIEW  ANGLE  BY  BAND  (PENDLETON,  CA) 
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Figure  76.  Histogram  Shift  with  View  Angle:  Pendleton,  CA  (Band  1) 
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Pendleton,  CA  -  Bond  2 


Figure  77.  Histogram  Shift  with  View  Angle:  Pendleton,  CA  (Band  2) 
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Pendleton,  CA  -  Bond  3 


Figure  78.  Histogram  Shift  with  View  Angle:  Pendleton,  CA  (Band  3) 
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Figure  79.  Histogram  Shift  with  View  Angle:  Pendleton,  CA  (Band  4) 
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Pendleton,  CA  -  Bond  5 


Figure  80.  Histogram  Shift  with  View  Angle:  Pendleton,  CA  (Band  5) 
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Figure  8 1 .  Histogram  Shift  with  View  Angle:  Pendleton,  CA  (Band  6) 
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Figure  82.  Histogram  Shift  with  View  Angle:  Pendleton,  CA  (Band  7) 


105 


Pendleton,  CA  -  Bond  8 


Figure  83.  Histogram  Shift  with  View  Angle:  Pendleton,  CA  (Band  8) 
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F. 


SPECTRAL  RADIANCE  DISPLAYING  VARIANCE  WITH  WAVELENGTH  AND  VIEW  ANGLE  (DUCK,  NC) 


Duck:  Ocean:  Spectral  Radiance 


Band 


-■-15:41:35 

-*-15:41:46 

-■-15:41:57 

-*-15:42:08 

-•-15:42:18 

-■-15:42:29 

-•-15:42:40 

-■-15:42:51 

-♦-15:43:01 

-*-15:43:12 

♦15:43:22 

-*-15:43:33 

15:43:43 

15:43:53 

15:44:03 


Figure  84.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Duck,  NC  (Ocean  ROI) 
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Duck:  Pond:  Spectral  Radiance 
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Figure  85.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Duck,  NC  (Pond  ROI) 
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Figure  88.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Duck,  NC  (Forest  ROI) 
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Duck:  Forest:  Spectral  Radiance 
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Figure  89.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Duck,  NC  (Golf  Course  ROI) 
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Duck:  Golf  Course:  Spectral  Radiance 
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Duck:  Parking  Lot:  Spectral  Radiance 
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Figure  91.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Duck,  NC  (Parking  Lot  ROI) 
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Duck:  Parking  Lot  2:  Spectral  Radiance 
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Figure  92.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Duck,  NC  (Parking  Lot  2  ROI) 


115 


Duck:  Roof  Top:  Spectral  Radiance 
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Figure  93.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Duck,  NC  (Roof  Top  ROI) 
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Figure  94.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Duck,  NC  (Roof  Top  2  ROI) 
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Duck:  Roof  Top  3:  Spectral  Radiance 
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Figure  95.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Duck,  NC  (Roof  Top  3  ROI) 
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G. 


SPECTRAL  RADIANCE  DISPLAYING  VARIANCE  WITH  WAVELENGTH  AND  VIEW  ANGLE  (PENDLETON, 
CA) 


Pendleton:  Ocean:  Spectral  Radiance 
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Figure  96.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Pendleton,  CA  (Ocean  ROI) 
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Pendleton:  Harbor:  Spectral  Radiance 
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Figure  97.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Pendleton,  CA  (Harbor  ROI) 
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Pendleton:  Marsh  (Veg):  Spectral  Radiance 
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Figure  98.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Pendleton,  CA  (Marsh  (Veg)  ROI) 
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Pendleton:  Field  (Veg):  Spectral  Radiance 
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Figure  99. 


Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Pendleton,  CA  (Field  (Veg)  ROI) 
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Pendleton:  Soil  (Wet):  Spectral  Radiance 
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Figure  100.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Pendleton,  CA  (Soil  (Wet)  ROI) 
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Pendleton:  Soil  (Dry):  Spectral  Radiance 
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Figure  101.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Pendleton,  CA  (Soil  (Dry)  ROI) 
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Pendleton:  Sand:  Spectral  Radiance 
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Figure  102.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Pendleton,  CA  (Sand  ROI) 
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Pendleton:  Sand  2:  Spectral  Radiance 
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Figure  103.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Pendleton,  CA  (Sand  2  ROI) 
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Pendleton:  Sand  3:  Spectral  Radiance 
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Figure  104.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Pendleton,  CA  (Sand  3  ROI) 
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Pendleton:  Tarmac:  Spectral  Radiance 
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Figure  105.  Spectral  Radiance  Displaying  Variance  with  Wavelength  and  View  Angle:  Pendleton,  CA  (Tarmac  ROI) 
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H. 


BRDF  RESULTS  BY  SEPARATION  ANGLE  FOR  REGIONS  OF 
INTEREST  (DUCK,  NC) 
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Figure  106.  BRDF  Results  by  Separation  Angle:  Duck,  NC  (Ocean  ROI) 
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Duck:  Pond:  BRDF 
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Figure  107.  BRDF  Results  by  Separation  Angle:  Duck,  NC  (Pond  ROI) 
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BRDF  (1/sr) 
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Figure  108.  BRDF  Results  by  Separation  Angle:  Duck,  NC  (Soil  (Dry)  ROI) 
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Duck:  Soil  (Wet):  BRDF 
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Figure  109.  BRDF  Results  by  Separation  Angle:  Duck,  NC  (Soil  (Wet)  ROI) 
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Duck:  Forest  BRDF 
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Figure  110.  BRDF  Results  by  Separation  Angle:  Duck,  NC  (Forest  ROI) 
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Duck:  Golf  Course:  BRDF 
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Figure  111.  BRDF  Results  by  Separation  Angle:  Duck,  NC  (Golf  Course  ROI) 
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Duck:  Baseball  Field:  BRDF 
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Figure  1 12.  BRDF  Results  by  Separation  Angle:  Duck,  NC  (Baseball  Field  ROI) 
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Duck:  Parking  Lot:  BRDF 
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Figure  113.  BRDF  Results  by  Separation  Angle:  Duck,  NC  (Parking  Lot  ROI) 
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BRDF  (1/sr) 


Duck:  Parking  Lot  2:  BRDF 
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Figure  1 14.  BRDF  Results  by  Separation  Angle:  Duck,  NC  (Parking  Lot  2  ROI) 
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Duck:  Roof  Top  2:  BRDF 
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Figure  1 16.  BRDF  Results  by  Separation  Angle:  Duck,  NC  (Roof  Top  2  ROI) 
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I. 


BRDF  RESULTS  BY  SEPARATION  ANGLE  FOR  REGIONS  OF 
INTEREST  (PENDLETON,  CA) 
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Figure  1 1 8.  BRDF  Results  by  Separation  Angle:  Pendleton,  CA  (Ocean  ROI) 
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Pendleton:  Harbor:  BRDF 
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Figure  119.  BRDF  Results  by  Separation  Angle:  Pendleton,  CA  (Harbor  ROI) 
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Pendleton:  Marsh  (Veg):  BRDF 
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Figure  120.  BRDF  Results  by  Separation  Angle:  Pendleton,  CA  (Marsh  (Veg)  ROI) 
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Pendleton:  Field  (Veg):  BRDF 
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Figure  121.  BRDF  Results  by  Separation  Angle:  Pendleton,  CA  (Field  (Veg)  ROI) 
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Pendleton:  Soil  (Wet):  BRDF 
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Figure  122.  BRDF  Results  by  Separation  Angle:  Pendleton,  CA  (Soil  (Wet)  ROI) 
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Pendleton:  Soil  (Dry):  BRDF 
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Figure  123.  BRDF  Results  by  Separation  Angle:  Pendleton,  CA  (Soil  (Dry)  ROI) 
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Pendleton:  Sand:  BRDF 
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Figure  124.  BRDF  Results  by  Separation  Angle:  Pendleton,  CA  (Sand  ROI) 
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Pendleton:  Sand  3:  BRDF 

0.0G 


0.05 


0.04 

♦  Band  1 
-■—Band  2 
Band  3 
Band  4 
Band  5 
Band  6 
— Band  7 
- Band  8 

0.01 


0 

-90  -85  -80  -75  -70  -65  -60  -55  -50 

Separation  Angle 

Figure  126.  BRDF  Results  by  Separation  Angle:  Pendleton,  CA  (Sand  3  ROI) 
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Pendleton:  Tarmac:  BRDF 
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Figure  127.  BRDF  Results  by  Separation  Angle:  Pendleton,  CA  (Tarmac  ROI) 
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